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Abstract. Let M and N be connected manifolds without boundary with 
dim(M) < dim(Af), and let M compact. Then shape space in this work is 
either tlic manifold of submanifolds of that are diffeomorphic to M, or 
the orbifold of unparametrized immersions of M in N. We investigate the 
Sobolcv Riemannian metrics on shape space; These are induced by metrics of 
the following form on the space of immersions: 

Gf{h,k)= f g{pi'h,k)vol{f*g) 
Jm 

where g is some fixed metric on N, f*g is the induced metric on M, h,k& 
r{f*TN) axe tangent vectors at / to the space of embeddings or immersions, 
and is a positive, selfadjoint, bijective scalar pseudo differential operator of 
order 2p depending smoothly on /. We consider later specifically the operator 
= 1 + j4A'', where A is the Bochner-Lapla<;ian on M induced by the metric 
f*g. For these metrics we compute the geodesic equations both on the space 
of immersions and on shape space, and also the conserved momenta arising 
from the obvious symmetries. We also show that the geodesic equation is well- 
posed on spaces of immersions, and also on difleomorphism groups. We give 
examples of numerical solutions. 



1. Introduction 

Many procedures in science, engineering, and medicine produce data in the form 
of shapes. If one expects such a cloud to foUow roughly a submanifold of a certain 
type, then it is of utmost importance to describe the space of all possible submani- 
folds of this type (we call it a shape space hereafter) and equip it with a significant 
metric which is able to distinguish special features of the shapes. Most of the met- 
rics used today in data analysis and computer vision are of an ad-hoc and naive 
nature; one embeds shape space in some Hilbert space or Banach space and uses 
the distance therein. Shortest paths are then line segments, but they leave shape 
space quickly. 

Riemannian metrics on shape space itself are a better solution. They lead to 
geodesies, to curvature and diffusion. Eventually one also needs statistics on shape 
space like means of clustered subsets of data (called Karcher means on Riemannian 



Date: June 27, 2011. 

2000 Mathematics Subject Classification. Primary 58B20, 58D15, 58E12. 
All authors were supported by FWF Project 21030 and by the NSF-Focused Research Group: 
The geometry, mechanics, and statistics of the infinite dimensional shape manifolds. 



2 



MARTIN BAUER, PHILIPP HARMS, PETER W. MICHOR 



manifolds) and standard deviations. Here curvature will play an essential role; 
statistics on Riemannian manifolds seems hopelessly underdeveloped just now. 

1.1. The shape spaces used in this work. Thus, initially, by a shape we mean 
a smoothly embedded surface in N which is diffeomorphic to M. The space of 
these shapes will be denoted — Be {M, N) and viewed as the quotient (see [T5] 
for more details) 

Be{M,N) = Emb(M, A^)/Diff(M) 

of the open subset Emb(Af, A^) C C°°{M,N) of smooth embeddings of M in N, 
modulo the group of smooth difFeomorphisms of M. It is natural to consider all 
possible immersions as well as embeddings, and thus introduce the larger space 
Bi = Bi{M,N) as the quotient of the space of smooth immersions by the group 
of difFeomorphisms of M (which is, however, no longer a manifold, but an orbifold 
with finite isotropy groups, see [T5]'). 

Emb(M, N) ^ Emb(M, N)/DiS{M) Be{M, N) 

P P P 

III 
Imm(Af, N) ^ Imm(M, N) / Diff (Af ) B,{M, N) 

More generally, a shape will be an element of the Cauchy completion (i.e., the 
metric completion for the geodesic distance) of Bi{M, N) with respect to a suitably 
chosen Riemannian metric. This will allow for corners. In practice, discretization 
for numerical algorithms will hide the need to go to the Cauchy completion. 

1.2. Where this work comes from. In [201, Michor and Mumford have investi- 
gated a variety of Riemannian metrics on the shape space 

Bi{S\R^) = lmm{S^,E.'^)/mS{S^) 

of unparametrized immersion of the circle into the plane. In |19l section 3.10] they 
found that the simplest such metric has vanishing geodesic distance; this is the 
metric induced by L^(arc length) on Imm(S'^, M'^): 

G°{h,k) = J^^{h{e),k{e)}\fmde, 

f e Imm(S'\M2)^ /j^fc g C°°{S\R'^) = T/Imm(S'\ E^). 

In [TS] they found that the vanishing geodesic distance phenomenon for the L^- 
metric occurs also in the more general shape space Imm(Af , iV)/ Diff (Af) where 
is replaced by a compact manifold M and Euclidean is replaced by Riemannian 
manifold N; it also occurs on the full diffeomorphism group Diff (iV), but not on the 
subgroup Diff (A'^j vol) of volume preserving diffeomorphisms, where the geodesic 
equation for the L^-metric is the Euler equation of an incompressible fluid. In 
pPi sections 3, 4 and 5] three classes of metrics were investigated: Almost local 
metrics on planar curves, Sobolev metrics on planar curves, and metrics induced 
from Sobolev metrics on the diffeomorphism group of the plane. The results about 
almost local metrics from [501 section 3] were generalized by the authors to the case 
of surfaces in [2]. 
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Now we take up the investigations from [20l section 4]. The immersion- Sobolev 
metric considered there is 

G^~ih,k)^ f Uh,k)+A.{Dlh,Dlk)).ds 
= / {Lp{h),k)ds where 
Lp{h) or Lpj{h) = {l + {-l)PA.D^/){h) and = 

\J9\ 

The interesting special case p — I and ^ — > oo has been studied in [551 US] and in 
|24) where an isometry to an infinite dimensional Grassmannian with the Fubini- 
Study metric was described. In this case, the metric reduces to: 

G)"''^-'-°^(h,k)^ f {D,{h),DM).ds 

The cases p = 1,2 and A — > oo have also been treated in [15], where estimates on 
the geodesic distance are proven and the metric completion of the space of curves 
is characterized. 

In this work we generalize the immersion-Sobolev metrics from |201 section 4] 
to higher dimensions and to non-flat ambient space, namely to the shape space 
Bi{M, N) = Imm(M, N)/DiS{M) of surfaces of type M in A^; here M is a compact 
orientable connected manifold of smaller dimension than N, for example a sphere 
S"',m< dim(iV). 

1.3. Riemannian metrics. The tangent space Tflmm{M, N) at an immersion / 
consists of all vector fields along /: 

T/Imm(M, N) = r{f*TN) ^ {h e C°^{M, TN) : tttn oh^ f}. 

A Riemannian metric on Imm(A/, N) is a family of positive definite inner products 
Gf{h, k) where / e Imm(A/, N) and h,k G T/Imm(M, N). Each metric is weak in 
the sense that G/, viewed as linear map from T/Imm(M, N) into its dual consisting 
of distributional sections of f*TN is injective. (But it can never be surjective.) 
We require that our metrics will be invariant under the action of Diff(M), hence 
the quotient map dividing by this action will be a Riemannian submersion. This 
means that the tangent map of the quotient map Imm(Af, A^) Bi{M,N) is a 
metric quotient mapping between all tangent spaces. Thus we will get Riemannian 
metrics on Bi. For any / e Imm{M,N) those vectors in Tflmm(M, N) which are 
G/-perpendicuar to the Diff(Af)-orbit through / are called horizontal (with respect 
to G). They form the G/-orthogonal space to the orbit. A priori we do not know 
that it is a complementary space. For the metrics considered in this work it will 
turn out to be a complement. 

The simplest inner product on the tangent bundle to Imm(Af , N) is 

G°(/i,fc)- / g{h,k)vo\{r9), 
Jm 

where g = { , ) is the Euclidean inner product on N. Since the volume form 
vo\{f*g) reacts equivariantly to the action of the group Diff(Af), this metric is 
invariant, and the map to the quotient B^ is a Riemannian submersion for this 
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metric. The G^-horizontal vectors in Tflmm{M, N) are just those vector fields 
along / which are pointwise ^-normal to f{M); we will call them normal fields. 

All of the metrics we will look at will be of the form (see section |6| : 

G^{h,k) = / g{pfh,k)vol{rg) 

where : T^Imm — > Tflmm is a positive bijective operator depending smoothly 
on /, which is selfadjoint unbounded in the Hilbert space Tyimm with inner product 
. We will assume that P is in addition equivariant with respect to reparanietriza- 
tions, i.e. 

pf°v opf o [ip-^y = Lp*{pf) for aU LP e Diff(M). 

The G^-horizontal vectors will be those h e T/Emb(Af,iV) = C°°{M,N) such 
that P^h is normal. 

The tangent map of the quotient map Emb(Af, N) — > Bi{M, N) is then an isom- 
etry when restricted to the horizontal spaces, just as in the finite dimensional situ- 
ation. Riemannian submersions have a very nice effect on geodesies: the geodesies 
on the quotient space Bi are exactly the images of the horizontal geodesies on the 
top space Imm; by a horizontal geodesic we mean a geodesic whose tangent lies in 
the horizontal bundle. The induced metric is invariant under the action of Diff (Af) 
and therefore induces a unique metric on Bi. See for example [2J section 1]. Later 
in section [s] we shall consider the special case P^ = 1 + AA^. 

1.4. Inner versus outer metrics. The metrics studied in this work are induced 
from Imm(M, N) on shape space. One might call them inner metrics since the 
differential operator governing the metric is defined intrinsic to M. Intuitively, 
these metrics can be seen as describing some elastic or viscous behaviour of the 
shape itself. 

In contrast to these metrics, there are also metrics induced from Diff(A^) on shape 
space. (The widely used LDDMM algorithm uses such a metric.) The differential 
operator governing these metrics is defined on all of N, even outside of the shape. 
Intuitively, these metrics can be seen as describing some elastic or viscous behaviour 
of the ambient space N that gets deformed as the shape changes. One might call 
these metrics outer metrics. 

1.5. Contributions of this work. 

• This work is the first to treat Sobolev inner metrics on spaces of immersed 
surfaces and on higher dimensional shape spaces. 

• It contains the first description of how the geodesic equation can be for- 
mulated in terms of gradients of the metric with respect to itself when the 
ambient space is not fiat. To achieve this, a covariant derivative on some 
bundles over immersions is defined. This covariant derivative is induced 
from the Levi-Civita covariant derivative on ambient space. 

• The geodesic equation is formulated in terms of this covariant derivative. 
Well-posedness of the geodesic equation is shown under some regularity as- 
sumptions that are verified for Sobolev metrics. Well-posedness also follows 
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for the geodesic equation on diffeomorphism groups, where this result has 
not yet been obtained in that full generality. 

• To derive the geodesic equation, a variational formula for the Laplacian 
operator is developed. The variation is taken with respect to the metric on 
the manifold where the Laplacian is defined. This metric in turn depends 
on the immersion inducing it. 

• It is shown that Sobolev inner metrics separate points in shape space when 
the order of the differential operator governing the metric is high enough. 
(The metric needs to be as least as strong as the iJ^-metric.) Thus Sobolev 
inner metrics overcome the degeneracy of the i^-metric. 

• The path-length distance of Sobolev inner metrics is compared to the 
Frechet distance. It would be desirable to bound Fechet distance by some 
Sobolev distance. This however remains an open problem. 

• Finally it is demonstrated in some examples that the geodesic equation for 
the 77^-metric on shape space of surfaces in M.^ can be solved numerically. 

Big parts of this work can also be found, partly in more details, in the doctoral 
theses of Martin Bauer ^ and Philipp Harms [U]. 

2. Content of this work 

This work progresses from a very general setting to a specific one in three steps. 
In the beginning, a framework for general inner metrics is developed. Then the 
general concepts carry over to more and more specific inner metrics. 

• First, shape space is endowed with a general inner metric, i.e with a metric 
that is induced from a metric on the space of immersions, but that is 
unspecified otherwise. It is shown how various versions of the geodesic 
equation can be expressed using gradients of the metric with respect to itself 
and how conserved quantities arise from symmetries. (This is section |4]) 

• Then it is assumed that the inner metric is defined via an elliptic pseudo- 
differential operator. Such a metric will be called a Sobolev-type metric. 
The geodesic equation is formulated in terms of the operator, and existence 
of horizontal paths of immersions within each equivalence class of paths is 
proven. (This is section |6]) Then estimates on the path-length distance 
are derived. Most importantly it is shown that when the operator involves 
high enough powers of the Laplacian, then the metric does not have the 
degeneracy of the L^-metric. (This is section [?]) 

• Motivated by the previous results it is assumed that the elliptic pseudo- 
differential operator is given by the Laplacian and powers of it. Again, the 
geodesic equation is derived. The formulas that are obtained are ready to 
be implemented numerically. (This is section [sj) 

The remaining sections cover the following material: 

• Section |3] treats some differential geometry of surfaces that is needed in 
this work. It is also a good reference for the notation that is used. The 
biggest emphasis is on a rigorous treatment of the covariant derivative. 
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Some material like the adjoint covariant derivative is not found in standard 
text books. 

Section [5] contains formulas for the variation of the metric, volume form, 
covariant derivative and Laplacian with respect to the immersion inducing 
them. These formulas are used extensively later. 

Section |9] covers the special case of flat ambient space. The geodesic equa- 
tion is simplified and conserved momenta for the Euclidean motion group 
are calculated. Sobolev-type metrics are compared to the Frechet metric 
which is available in flat ambient space. 

Section 10 treats diffeomorphism groups of compact manifolds as a special 
case of the theory that has been developed so far. 

In section [TT] it is shown in some examples that the geodesic equation on 
shape space can be solved numerically. 



3. Differential geometry of surfaces and notation 

In this section the differential geometric tools that are needed to deal with im- 
mersed surfaces are presented and developed. The most important point is a rig- 
orous treatment of the covariant derivative and related concepts. 

The notation of 17J is used. Some of the definitions can also be found in [T^]. 
A similar exposition in the same notation is ^ . 

3.1. Basic assumptions and conventions. 

Assumption. It is always assumed that M and N are connected manifolds of finite 
dimensions m and n, respectively. Furthermore it is assumed that M is compact, 
and that N is endowed with a Riemannian metric g. 

In this work, immersions of A/ into N will be treated, i.e. smooth functions M — > 
N with injective tangent mapping at every point. The set of all such immersions 
will be denoted by Imm(M, N). It is clear that only the case dim(il/) < dim(A^) is 
of interest since otherwise Imm(M, N) would be empty. 

Immersions or paths of immersions are usually denoted by /. Vector fields on 
Imm(M, N) or tangent vectors with foot point /, i.e., vector fields along /, will be 
called h, k, m, for example. Subscripts like ft — dtf = df /dt denote differentiation 
with respect to the indicated variable, but subscripts are also used to indicate the 
foot point of a tensor field. 

3.2. Tensor bundles and tensor fields. The tensor bundles 

TIM TlM®f*TN 

Y y 

M M 

will be used. Here TJM denotes the bundle of ( ^ )-tensors on Af, i.e. 

r s 

TIM = (^TM®<^ T*M, 
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and f*TN is the pullback of the bundle TN via /, see [TT!, section 17.5]. A tensor 
field is a section of a tensor bundle. Generally, when E is a bundle, the space of its 
sections will be denoted by r(£'). 

To clarify the notation that will be used later, some examples of tensor bun- 
dles and tensor fields are given now. S^T* M = L^y^{TM;R) and K^T* M = 
L^u(TM;R) are the bundles of symmetric and alternating (°)-tensors, respec- 
tively. f]'=(M) = r(A'=T*M) is the space of differential forms, X(M) = r{TM) is 
the space of vector fields, and 

r{f*TN) = {he C°°{M, TN) : TTAT o /i = /} 

is the space of vector fields along f. 

3.3. Metric on tensor spaces. Let g e r(S'?.gT*7V) denote a fixed Riemannian 
metric on N. The metric induced on M by f € Imm(M, N) is the pullback metric 

g^rge T{Si,T*M), g{X, Y) - Y) = 9{Tf.X, TJ.Y), 

where X, Y are vector fields on M . The dependence of g on the immersion / should 
be kept in mind. Let 

b = 5 : TM T*M and jj = : r*M ^ TM. 

g can be extended to the cotangent bundle T*M = T^M by setting 

g-\a,l3)=gUa,l3) = a{l3^) 

for a, /3 e T*M. The product metric 

extends g to all tensor spaces T^M, and gl^g yields a metric on TJM f*TN. 

3.4. Traces. The trace contracts pairs of vectors and co-vectors in a tensor prod- 
uct: 

Tr : T*M (g) TM = L{TM, TM) ^ M x M 

A special case of this is the operator ix inserting a vector X into a co-vector or 
into a covariant factor of a tensor product. The inverse of the metric g can be used 
to define a trace 

Tr» : T*M ® T* M ^ M xR 

contracting pairs of co-vecors. Note that Tr^ depends on the metric whereas Tr 
does not. The following lemma will be useful in many calculations: 

Lemma. 

gO(B, C) = TT{g-^Bg-^C) for B,C e T^M if B or C is symmetric. 
(In the expression under the trace, B and C are seen as maps TM — > T* M .) 
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Proof. Express everything in a local coordinate system u^, . . . , m™ of M. 

ik jl 

^Y.9''B.k9'''Cji =Y,g^'B,k9'''Ci, = Tv{g-' Bg-^C) 

ijkl ijkl 

Note that only the symmetry of C has been used. □ 

3.5. Volume density. Let Vol(M) be the density bundle over M, see [T71 sec- 
tion 10.2]. The volume density on AI induced by / G Imm(M, iV) is 

vol(g)=vol(rg)er(Vol(M)). 

The volume of the immersion is given by 

Vol(/) = / vol(/*g) = / vol(.g). 
Jai Jm 

The integral is well-defined since M is compact. If M is oriented the volume density 
may be identified with a differential form. 

3.6. Metric on tensor fields. A metric on a space of tensor fields is defined by 
integrating the appropriate metric on the tensor space with respect to the volume 
density: 

gl{B,C)= f g:{Bix),CXx))Yo\ig)ix) 

for B,C e r(r;M), and 

^g{B,C)= f g:®g{B{x),C{x))yo\{g){x) 

for B,C e r(TJM(g f*TN), f e Imm(Af,iV). The integrals are well-defined 
because M is compact. 

3.7. Covariant derivative. Covariant derivatives on vector bundles as explained 
in [17l sections 19.12, 22.9] will be used. Let V^, be the Levi-Civita covariant 
derivatives on {M,g) and {N,g), respectively. For any manifold Q and vector field 
X on Q, one has 

: C°°(Q,rAf) -> C°°(Q,rM), h^V^^h 

V\:C°°{Q,TN)^C°°{Q,TN), h^V\h. 

Usually the symbol V will be used for all covariant derivatives. It should be kept in 
mind that depends on the metric g ~ f*g and therefore also on the immersion 
/. The following properties hold [T71 section 22.9]: 

(1) ^ X respects base points, i.e. tt o Vx^ — t: oh, where tt is the projection of 
the tangent space onto the base manifold. 

(2) y xh is C°°-linear in X. So for a tangent vector X^ G T^Q, ^x^h makes 
sense and equals {\/xh){x)- 

(3) Vxh is M-linear in h. 

(4) Vxia.h) = da{X).h -t- a.V xh for a £ C°°{Q), the derivation property of 
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(5) For any manifold Q and smooth mapping q : Q ^ Q and Yy e TyQ one 
has ^Tq.Yyh = ^Yy{hoq). If F e ^{Qi) and X e are g-related, then 
Vy(/i o q) = (Vxh) o q. 

The two covariant derivatives and can be combined to yield a covariant 
derivative Vx acting on C°°{Q, M (^TN) by additionally requiring the following 
properties JTi section 22.12]: 

(6) Vx respects the spaces C°°{Q,T^M ® TN). 

(7) y x{h ® k) = {V xh) (E) k + h iE> (Vjcfc), a derivation with respect to the 
tensor product. 

(8) Vx commutes with any kind of contraction (see \T7l section 8.18]). A 
special case of this is 

\/x{a(Y)) = i^xaW) + oti^ x^) for a (g) F : iV -> TIM. 

Property ([T]) is important because it implies that Vx respects spaces of sections of 
bundles. For example, for Q — M and / e C°°{M, N), one gets 

Vx : r(T;M ® f*TN) r{T^M «) f*TN). 

3.8. Swapping covariant derivatives. Some formulas allowing to swap covariant 
derivatives will be used repeatedly. Let / be an immersion, h a vector field along 
/ and X, Y vector fields on M. Since V is torsion-free, one has TT, section 22.10]: 

(1) VxTf.Y - VyTf.X - Tf. [X, Y] - Tor(T/.X, Tf.Y) = 0. 
Furthermore one has |T71 section 24.5]: 

(2) VxVyh - Vy^xh - Vix,Y]h = W> o (T/.X, Tf.Y)h, 
where EF G Vl^ (^N ; L{T N ,T N)) is the Riemann curvature tensor of {N,g). 

These formulas also hold when / :IRxM— >A^isa path of immersions, h : 
R X M — )■ TN is a vector field along / and the vector fields are vector fields on 
M X Af. A case of special importance is when one of the vector fields is (9t, Om) and 
the other (OR,y), where y is a vector field on M. Since the Lie bracket of these 
vector fields vanishes, ([T]) and ^ yield 

(3) V(a,^OM)T^/-(0K,>^) - V(o.,y)T/.(a„OM) = 
and 

(4) V(a,,OM)V(o.,y)/^- V(o«,y)V(a„o.,,)/i = i?^(T/.(at, Om), r/.(OR, 

3.9. Second and higher covariant derivatives. When the covariant derivative 
is seen as a mapping 

V : r(r;Af ) t{ti^^m) or V : r(T;Af ® f*TN) r(r;+iA/ «) f*TN), 

then the second covariant derivative is simply VV — V'^. Since the covariant 
derivative commutes with contractions, can be expressed as 

Vx,y — i^YixV^ = tyVxV = VxVy - Vy^y for X,Y £ X{M). 

Higher covariant derivates are defined accordingly as V*^, k > 0. 
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3.10. Adjoint of the covariant derivative. The covariant derivative 

V : t{t;m) ^ r(T;+iM) 

admits an adjoint 

V* : r(r;+iA/) r(T;M) 

with respect to the metric g, i.e.: 

In the same way, V* can be defined when V is acting on r(TJM (g) f*TN). In 
either case it is given by 

V*B = -Tr^'(VB), 

where the trace is contracting the first two tensor slots of WB. This formula will 
be proven now: 

Proof. The result holds for decomposable tensor fields (3 ^ B E T{T^^iM) since 
g-(v*{P(^B),c) ^g^^((3®B,Vc) ^gl(B,Vp,c) 

= [ Cp,g:iB,C)vo\ig)^ [ g^iV p,B,C)vo\{g) 

JM JM 

- / -g:{B,C)C,3,vo\{g)- f g:( Tr^ (/? ® VB), C)vo%) 

JM JM 

= - div(/3«)B - Ti^ip (g) WB),C^ 

= 5j(-div(/3«)B + Tr^'((V^)(g)B) -Tr^'(V(/3(8)S)),c) 

= - div(/3«)B + Tr^'(V/3)B - Tr^ (V(/3 ® B)), c) 

= .gr(0-TrS(V(/3®S)),c) 

Here it has been used that = 0, that commutes with any kind of con- 

traction and acts as a derivation on tensor products 17, section 22.12] and that 
div(X) = Tr(VX) for all vector fields X [T71 section 25.12]. To prove the result for 
13® B € r(rj+iM (g) f*TN) one simply has to replace gl by gl^g. □ 

3.11. Laplacian. The definition of the Laplacian used in this work is the Bochner- 
Laplacian. It can act on all tensor fields B and is defined as 

AB = y*VB = - Tr^iV^B). 

3.12. Normal bundle. The normal bundle Nor(/) of an immersion / is a sub- 
bundle of f*TN whose fibers consist of all vectors that are orthogonal to the image 
of/: 

Nor(/), = {r e T^(,)iV : VX e r,M : g{Y,TJ.X) = O}. 

If dim(Af) = dim(Af) then the fibers of the normal bundle are but the zero vector. 
Any vector field h along / G Imm can be decomposed uniquely into parts tangential 
and normal to / as 

h = Tf.h^ + h^, 

where is a vector field on M and is a section of the normal bundle Nor(/). 
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3.13. Second fundamental form and Weingarten mapping. Let X and Y be 

vector fields on M. Then the covariant derivative WxTf.Y splits into tangential 
and a normal parts as 

VxTf.Y = Tf.{VxTf.YY + (VxTf.Y)^ = Tf.VxY + S{X,Y). 

S is the second fundamental form of f. It is a symmetric bilinear form with values 
in the normal bundle of /. When Tf is seen as a section of T*M (g) f *TN one has 
S = VTf since 

S{X, Y) = VxTf.Y - Tf.VxY = {VTf){X, Y). 
The trace of S is the vector valued mean curvature Tr^(S') G r(Nor(/)). 

4. Shape space 

Briefly said, in this work the word shape means an unparametrized surface. (The 
term surface is used regardless of whether it has dimension two or not.) This sec- 
tion is about the infinite dimensional space of all shapes. First some spaces of 
parametrized and unparametrized surfaces are described, and it is shown how to 
define Riemannian metrics on them. The geodesic equation and conserved quanti- 
ties arising from symmetries are derived. 

The agenda that is set out in this section will be pursued in section [6] when the 
arbitrary metric is replaced by a Sobolev-type metric involving a pseudo-differential 
operator and later in section [8] when the pseudo-differential operator is replaced by 
an operator involving powers of the Laplacian. 

4.1. Riemannian metrics on immersions. The space of smooth immersions of 
the manifold M into the manifold N will be denoted by Imm(M, N) or briefly Imm. 
It is a smooth Frechet manifold containing the space Emb(M, N) of embeddings 
of M into N as an open subset [111 theorem 44.1]. Consider the following natural 
bundles of k-multilinear mappings: 

L'= (Tlmm; M) L'' (TImm; Tlmm) 

Imm Imm 
These bundles are isomorphic to the bundles 

L ( ^'^Tlmm; M ) L ( ^'^Tlmm; Tlmm ) 



Imm Imm 

where denotes the c°°-completed bornological tensor product of locally convex 
vector spaces [HI section 5.7, section 4.29]. Note that L(TImm;TImm) is not 
isomorphic to T*Imm <§ Tlmm since the latter bundle corresponds to multilinear 
mappings with finite rank. 
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It is worth to write down more explicitly what some of these bundles of multi- 
linear mappings are. The tangent space to Imm is given by 

T/Imm = Cf{M,TN) {h G C°°{M,TN) : Tr^v o = /}, 

TImm = Cf^^„^{M,TN) := {h e C°°{M,TN) : ttn o h e Imm}. 

Thus Tflmm is the space of vector fields along the immersion /. Now the cotangent 
space to Imm will be described. The symbol ®c°=(j\/) means that the tensor product 
is taken over the algebra C°°(M). 

r;imm = L(r/Imm;M) = Cf{M,TN)' = C°°(M)' ®c°-(M)Cf(M,T* N) 

r*Imm = L(rimm;M) = C°°(M)' ®c^(M)CZra{MX N) 

The bundle L^^j„(rimm; M) is of interest for the definition of a Riemannian metric 
on Imm. (The subscripts sym and ait indicate symmetric and alternating multilinear 
maps, respectively.) Letting ®s denotes the symmetric tensor product and ®s the 
c°°-completed bornological symmetric tensor product, one has 

Lly^{Tfliam-R) = (T^Imm is T/Imm)' ^ {Cf(M,TN) §)s Cf{M,TN))' 

= {Cf{M,TN ®s TN))' 

= C°°iMY ^c^^M)Cf{M,T*N T*N) 

L2^^(TImm;R) = C°°(M)' ^c^^M)CS^miM,T* N (^s T*N) 

A Riemannian metric G on Imm is a section of the bundle Lgy,jj(TImm; M) such 
that at every / € Imm, G/ is a symmetric positive definite bilinear mapping 

Gf : T/Imm x T/Imm M. 

Each metric is weak in the sense that G / , seen as a mapping 

Gf : Tjimm — > T^Imm 

is injective. (But it can never be surjcctive.) 

4.2. Covariant derivative on immersions. The covariant derivative de- 
fined in section [3.71 induces a covariant derivative over immersions as follows. Let 
Q be a smooth manifold. Then one identifies 

heG°° (g, TImm(Af, N)) and X e X(Q) 

with 

/i^ G C°°(Q X M, TN) and {X, Om) e X{Q x M). 

As described in section [3771 one has the covariant derivative 

Vf^^.^^^h^ eG^{QxM,TN). 

Thus one can define 

Vx/i= (vf^,o..)/i^)'' e G°^{Q,Tlmm{M,N)). 

This covariant derivative is torsion-free by section [sTsj formula Q. It respects the 
metric g but in general does not respect G. 

It is helpful to point out some special cases of how this construction can be 
used. The case Q = M will be important to formulate the geodesic equation. The 
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expression that will be of interest in the formulation of the geodesic equation is 
Vg^/t, which is well-defined when / : M — > Imm is a path of immersions and 
/( : M Tlmm is its velocity. 

Another case of interest is Q = Imm. Let h,k,m € X(Imm). Then the covariant 
derivative Vmh is well-defined and tensorial in m. Requiring Vm to respect the 
grading of the spaces of multilinear maps, to act as a derivation on products and to 
commute with compositions of multilinear maps, one obtains as in section [377| a co- 
variant derivative Vm acting on all mappings into the natural bundles of multilinear 
mappings over Imm. In particular, V mP and V^G are well-defined for 

P e r(L(TImm;rimm)), G G r(L^y„,(rimm; M)) 

by the usual formulas 

{V,nP){h) = Vm{P{h)) - P{Vmh), 

{W^G){h,k) = W,n{G{h,k)) - G{W„Ak) - G{h,V„,k). 

4.3. Metric gradients. The metric gradients H,K G r(L^ (Tlmm; Tlmm)) are 
uniquely defined by the equation 

(V,„G)(/i,fc) ^ G{K{h,m),k) = G{m,H{h,k)), 

where h,k,m are vector fields on Imm and the covariant derivative of the metric 
tensor G is defined as in the previous section. (This is a generalization of the 
definition used in [201 that allows for a curved ambient space N ^ M".) 

Existence of if, K has to proven case by case for each metric G, usually by partial 
integration. For Sobolev metrics, this will be proven in sections |8.2| and [8.3[ 

Assumption. Nevertheless it will be assumed for now that the metric gradients 
H, K exist. 

4.4. Geodesic equation on immersions. 

Theorem. Given H, K as defined in the previous section and V as defined in 
section \4-^ the geodesic equation reads as 

^dJt = \Hf{fuft)^Kf{fuft). 

This is the same result as in [201 section 2.4], but in a more general setting. 

Proof. Let / : (— e, e) x [0, 1] x M — TV be a one-parameter family of curves of 
immersions with fixed endpoints. The variational parameter will be denoted by 
s S (— e,e) and the time-parameter by i G [0, 1]. In the following calculation, let 
Gf denote G composed with /, i.e. 

G/ : M ^ Imm L,2y„^(TImm; M). 

Remember that the covariant derivative on Imm that has been introduced in sec- 
tion |4j2] is torsion- free so that one has 



^ajs - ^ajt = Tf.[dud,] + Tor(/t, /,) = 0. 
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Thus the first variation of the energy of the curves is 
1 /-i „ I '■^ 



ds^ Gf{ftJt)dt^- (Vg^Gf){fuft)+ I GfiVoJuft)dt 







\^f.G){ftJt)+ C Gf{VoJsJt)dt 
Jo 

{VfG){ft,ft)dt+ [ dt Gf{f,Jt)dt 
Jo 

\wf,G){fsJt)dt- r Gf{fs,W9jt)dt 
Jo 

- ^ g(/., ^H{fuft)+0~ Kifuft) - Vajt) 
If /(O, •, •) is energy-minimizing, then one has at s = that 



dt. 



^iJ(/„/0-^(/t,/t)-VaJt-0. □ 

4.5. Geodesic equation on immersions in terms of the momentum. In the 

previous section the geodesic equation for the velocity ft has been derived. In many 
applications it is more convenient to formulate the geodesic equation as an equation 
for the momentum G(/f , •) G T^Imm. G{ft, •) is an element of the smooth cotangent 
bundle, also called smooth dual, which is given by 

G(TImm) ]J G/(T/Imm) = ]J {Gf{h, ■) : h e T/Imm} C T*Imm. 

/tElmm /(Elnim 

It is strictly smaller than T'*Imm since at every / G Imm the metric G/ : Tyimm — > 
Tjlmm is injective but not surjective. It is called smooth since it does not contain 
distributional sections of f*TN, whereas T^Imm does. 

Theorem. The geodesic equation for the momentum p G T*Imm is given by 

P = Gift, •) 
\/o,P=^Gf{H{ft,ft),-), 



where H is the metric gradient defined in section \4.3\ and V is the covariant deriv- 
ative action on mappings into r*Imm as defined in section\4- 



Gf.R^ Imm L^^„(TImm; M). 



Proof. Let G / denote G composed with the path / : M Imm, i.e. 
n r ■ m ^ Tmm ^ r 

syn 



Then one has 

Va,p = Va, {Gfift, •)) - (Va,G/)(/„ •) + Gy(VaJt, •) 
= (V/,G)(/t, •) + Gf[^H{ft, ft) - K{fu ft), .) 

= Gf{K{ft,ft), ■) + Gf{^-H{ft,ft) - K{ft,ft), •) 



□ 
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This equation is equivalent to Hamilton's equation restricted to the smooth 
cotangent bundle: 

P = G{fu •) 
pt = (grad" E){p). 

Here uj denotes the restriction of the canonical symplectic form on r*Imm to the 
smooth cotangent bundle and E is the Hamiltonian 

E : G(TImm) ^ R, E{p) = G^^ip.p) 

which is only defined on the smooth cotangent bundle. 

4.6. Shape space. Diff(M) acts smoothly on Imm(M,7V) and Emb(M,iV) by 
composition from the right. For Imm, the action is given by the mapping 

Imm(M, N) x Diff(A/) ^ Imm(Af, N), (/, ip) ^ r{f, tp) = r'^(/) = / o tp. 

The tangent prolongation of this group action is given by the mapping 

rimm(Af, N) x Diff(Af ) ^ TImm(M, N), (h, tp) ^ Tr'^{h) = hotp. 

Shape space is defined as the orbit space with respect to this action. That means 
that in shape space, two mappings differing only in their parametrization will be 
regarded the same. 

Theorem. Let M be compact and of dimension < n. Then Emb(M, A'^) is the total 
space of a smooth principal fiber bundle with structure group Diff (M), whose base 
manifold is a Hausdorff smooth Frechet manifold denoted by 

BeiM, N) = Emb(A/, N)/DiS{M). 

However, the space 

Bi {M, N) = Imm(Af, N) /Diff (Af) 

is not a smooth manifold, but has singularities of orbifold type: Locally, it looks like 
a finite dimensional orbifold times an infinite dimensional Frechet space. 

The proof for immersions can be found in [B] and the one for embeddings in [TU 
section 44.1]. As with immersions and embeddings, the notation Bi, B^ will be used 
when it is clear that M and N are the domain and target of the mappings. 

4.7. Riemannian metrics on shape space. We start with a metric G on Imm. 
The mapping tt : Imm Bi is a submersion of smooth manifolds, that is, Ttt : 
Tlmm — > TBi is surjective. 

V = V{tt) := ker(r7r) C Tlmm 

is called the vertical subbundle. The horizontal subbundle is the G-orthogonal sub- 
space of V: 

Hor Hor(7r, G) := ^(Tr)-^ C Tlmm. 

It need not be a complement to V (recall that the metric is weak; the complement 
could be in a suitable completion of the tangent space). For all metrics in this 
paper it will turn out to be a complement, however. Then any vector h G Tlmm 
can be decomposed uniquely in vertical and horizontal components as 
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This definition extends to the cotangent bundle as follows: An element of r*Imm 
is called horizontal when it annihilates all vertical vectors, and vertical when it 
annihilates all horizontal vectors. 

In the setting described so far, the mapping 

TfTTlnorf ■ Hor/ -> T^(f)Bi 

is an isomorphism of vector spaces for all / G Imm. This isomorphism will be used 
to describe the tangent space to Bi. If both Imm and Bi are Ricmannian manifolds 
and if this isomorphism is also an isometry for all / G Imm, then vr is called a 
Riemannian submersion. In that case, the metric G on Imm is Diff(Af)-invariant. 
This means that G = {r'^)*G for all Lp € Diff (M), where denotes the right action 
of Lp on Imm that was described in section |4.6[ This condition can be spelled out 
in more details using the definition of as follows: 

Gf{h,k) = {{r^yG){h,k) ^ Gr.if){Tr^{h),Tr^{k)) = G f.^{h o ip,k o ip). 

The following theorem establishes the converse statement: 

Theorem. Given a DiS (M) -invariant Riemannian metric on Imm, there is a 
unique Riemannian metric on the quotient space Bi such that the quotient map 
TT : Imm Bi is a Riemannian submersion. 

Proof. If the horizontal bundle Hor/ is a complement to V/ then T/tt : Hor/ — > 
TT^(^f^Bi is an isomorphism (off the orbifold singularities of Bi) and we can induce 
the metric on TT^(^f)Bi which is independent of the choice of / in the fiber over 7r(/) 
by the the Diff (M)-invariance of the metric. If it is not a complement one has to 
consider the metric quotient norm. See for example [19, section 3]. □ 

Assumption. R will always be assumed that a DiS{M) -invariant metric G on 
Imm(M, N) is given and that shape space Bi is endowed with the unique metric 
such that the quotient map is a Riemannian submersion. 

4.8. Riemannian submersions and geodesies. It follows from the general the- 
ory of Riemannian submersions that horizontal geodesies in the top space corre- 
spond nicely to geodesies in the quotient space: 

Theorem. Let c : [0, 1] — !■ Imm be a geodesic. 

(1) If c'{t) is horizontal at one t, then it is horizontal at all t. 

(2) If c'(t) is horizontal then tt o c is a geodesic in Bi. 

(3) If every curve in Bi can be lifted to a horizontal curve in Imm, then there 
is a one-to-one correspondence between curves in Bi and horizontal curves 
in Imm. This implies that instead of solving the geodesic equation on Bi 
one can equivalently solve the equation for horizontal geodesies in Imm. 

See [iTJ section 26] for the proof. 

4.9. Geodesic equation on shape space. Theorem |4 . 8| applied to the Riemann- 
ian submersion vr : Imm — > Bi yields: 



SOBOLEV METRICS ON SHAPE SPACE OF SURFACES 



17 



Theorem. Assuming that every curve in Bi can be lifted to a horizontal curve in 
Imm, the geodesic equation on shape space is equivalent to 



1}°' e Hor 



(1) 



where f is a horizontal curve in Imm, where H, K are the metric gradients defined 
in section 4-3 and where V is the covariant derivative defined in section\4.S\ 



1 \ hor 



This is a consequence of the Diff (Af)-invariance of the metric G and the conser- 
vation of the reaparametrization momentum. A general proof can be found in 
section 3.14]. 

It wiU be shown in section [6^ that curves in Bi can be hfted to horizontal curves 
in Imm for the very general class of Sobolev type metrics. Thus all assumptions 
and conclusions of the theorem hold. 



4.10. Geodesic equation on shape space in terms of the momentum. As 

in the previous section, theorem |4.8| will be applied to the Riemannian submersion 
TT : Imm — > Bi. But this time, the formulation of the geodesic equation in terms of 
the momentum will be used, see section [431 will be seen in section [6. 11[ this is 
the most convenient formulation of the geodesic equation for Sobolev-type metrics. 

Theorem. Assuming that every curve in Bi can be lifted to a horizontal curve in 
Imm, the geodesic equation on shape space is equivalent to the set of equations 

( p=G/(/t,-)eHorcT*Imm, 

j(Va.Pr'- = ^G/(iI(/„/0,-f°^ 

Here f is a curve in Imm, H is the metric gradient defined in section \4-3[ and 
V is the covariant derivative defined in section \4-S\ f is horizontal because p is 
horizontal. 



5. Variational formulas 

Recall that many operators like 

g = f*g, S = S^, vo\{g), V = V^, A = A^, ... 

implicitly depend on the immersion /. In this section their derivative with respect 
to / which is called their first variation will be calculated . These formulas will be 
used to calculate the metric gradients that are needed for the geodesic equation. 

This section is based on [2], see also [11]. Some but not all of the formulas were 
known before [5j [18] . More variational formulas can be found in O [231 [1] ■ 



5.1. Paths of immersions. All of the differential-geometric concepts introduced 
in section |3] can be recast for a path of immersions instead of a fixed immersion. 
This allows to study variations of immersions. So let / : K — >■ Imm(M, N) be a 
path of immersions. By convenient calculus [2], / can equivalently be seen as 
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/ : M X A/ — > such that f{t, ■) is an immersion for each t. The bundles over M 
can be replaced by bundles over M x M: 

pr*TJAf pr* T^M(8>f*TN Nor(/) 

Y 1 

M X Af M X Af Rx M 

Here pr2 denotes the projection pr2 : M x Af — > Af . The covariant derivative W zh is 
now defined for vector fields Z on M x Af and sections h of the above bundles. The 
vector fields (9t,0M) and (Om,X), where X is a vector field on Af, are of special 
importance. In later sections they will be identified with dt and X whenever this 
does not pose any problems. Let 

inst : Af -> R X Af , x H> {t, x). 

Then by property |5] from section |3.7| one has for vector fields X, V on Af 

VxTfit, ■).¥ = VxT{f o inst) o Y = VxTf o Tins* oY 

= VxTf o (Or, Y) o inSt = VTins, oA-T/ o (Or, Y) 

= (V(o,,x)T/o(OR,r))oinst. 

This shows that one can recover the static situation at t by using vector fields on 
M X Af with vanishing K-component and evaluating at t. 

5.2. Directional derivatives of functions. The following ways to denote direc- 
tional derivatives of functions will be used, in particular in infinite dimensions. 
Given a function F{x,y) for instance, 

D(x,h)F will be written as a shorthand for dt\oF{x + th,y). 

Here (x, h) in the subscript denotes the tangent vector with foot point x and di- 
rection h. If F takes values in some linear space, this linear space and its tangent 
space will be identified. 

5.3. Setting for first variations. In all of this chapter, let / be an immersion 
and /( e T^Imm a tangent vector to /. The reason for calling the tangent vector 
/( is that in calculations it will often be the derivative of a curve of immersions 
through /. Using the same symbol / for the fixed immersion and for the path of 
immersions through it, one has in fact that 

Dif,n)F^dtF{f{t)). 

5.4. Variation of equivariant tensor fields. Let the mapping 

F : Imm(Af, N) r{T^M) 

take values in some space of tensor fields over Af , or more generally in any natural 
bundle over M, see [13] . 

Lemma. If F is equivariant with respect to pullbacks by diffeomorphisms of M , 
i.e. 

F{f) = {^*F){f) = v*(F{{^-^rf)) 
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for all If e Diff(A/) and f G Imm(A/, A^), then the tangential variation of F is its 
Lie-derivative: 

DiLTfjr^F = dt\oF[foFlf) = dt\oF[{Flfrf) 
= d,\o{Flfy{F{f))^C^.r{Fif)). 

This allows us to calculate the tangential variation of the puUback metric and 
the volume density, for example. 

5.5. Variation of the metric. 

Lemma. The differential of the pullback metric 

Imm -> r{SloT*M), 
_ f ^9 = 1*9 

is given by 

D^fj^^g = 2Sym5(V/„r/) = ~2g{ft,S) + 2SymV(/7)^ 
= -2g{ft,S)+Cf.g. 

Here Sym denotes the symmetric part of the tensor field C of type ( 2 ) given by 
( Sym(C)) {X, Y) := \ {C{X, Y) + C{Y, X)) . 

Proof. Let / : M x Af iV be a path of immersions. Swapping covariant derivatives 
3.8 formula (|3| one gets 



as m section 



dt{giX,Y)) ^ dt{g{Tf.X,Tf.Y)) =giVa^Tf.X,Tf.Y) +g{Tf.X,V9,Tf.Y) 
^g(VxfuTf.Y)+g{Tf.X,VYft) = {2Sym g{V fuT f)){X,Y). 
Splitting ft into its normal and tangential part yields 

2SymgiV fuTf) ^ 2SymgiV f^ +VTf.f J ,Tf) 

- -2Symg(/,^, VT/) + 2 Sym.g(V/7, •) 
= ~25(/,^,5) + 2SymV(/7)^ 

Finally the relation 

F^U,Tf.fJ)9 = 2SymV(/7)^ = C^jg 
follows either from the equivariance of g with respect to puUbacks by diffeomor- 



phisms (see section 5.4) or directly from 

{Cxg){Y, Z) = Cx{g{Y, Z)) - g{CxY, Z) ~ g{Y,CxZ) 

= Vx {g{Y, Z)) - g{VxY - VyX, Z) - g{Y, V^Z - VzX) 

- giVyX, Z) + g(Y, VzX) = {VyX)\Z) + (VzXfiY) 

= (VyX^)(Z) + (VzX")(y) = 2Sym(V(X^))(r,Z). □ 
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5.6. Variation of the inverse of the metric. 

Lemma. The differential of the inverse of the pullback metric 

Imm ^ r{L{T*M,TM)), 



f ^ = {1*9)-' 

is given by 

DUM9-' = DuMifg)-' = 2g{f,^,g-'Sg-') + Cfr{g-') 



Proof. 

dtg-' = -9-\dtg)g-' = -g-\- 2g{ft ,S) + C^jg)g-^ 

= 2g{ft^,g^^Sy j^^f^ 



= 2g-^g{ft,S)g-' -g-\Cfjg)g-^ = 2g{f,^,g-^Sg-') + CfAg-^) □ 



5.7. Variation of the volume density. 
Lemma. The differential of the volume density 

Imm Vol(Af), 



/ ^ vol(g) = vol(/*g) 
is given by 

D^jj^jY0\{g) = Tr^ (g(V/„ T/)) vol(5) = ( div''(/7) - Tr^(5))) vol(5). 

Proof. Let g{t) G r{SyQT* M) be any curve of Riemannian metrics. Then 

dtyo\{g) = ^TT{g'\dtg)vo\{g). 
This follows from the formula for vol(5) in a local oriented chart {u^, . . . u™) on M: 
dtvo\{g) = dt^ det{{g^j)^j) du'^ A ■ ■ ■ A du^ 

= , , , Tr(adj(g)a,g) du' A ■ ■ ■ A 

2^det((g,j),y) 

^ Ti{Aei{{g,j)^j)g-^dtg) du^ A ■ ■ ■ A du"" 



2^det((g,j),j) 
\TT{g-\dtg)^o\{g) 



Now one can set g = f* g and plug in the formula 

dtg^dt{f*g)^2Symg{yft,Tf) 
from |5.5] This immediately proves the first formula: 

5tVol(5) - i Tr (.g-i.2Symg(V/t,r/)) = Tr^ {g{^ft,Tf)). 
Expanding this further yields the second formula: 
9*vol(.g) - Tr^ [^^gifuTf] - g{fu VT/)) 

= Tr^ {VgifuTf) - a{fu S)) - -V*g{fu Tf) - g{fuTT^{S)) 
= -V*{{fJf)~g{ft.Tr^{S))^Aiy{fJ)-g{ft,Tr^{S)). 
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Here it has been used that 

VT/ = S and div(/7) = Tr(V/7) = Tr^ ((V/7)^) = -V* {(f^f). 

Note that by |5.4[ the formula for the tangential variation would have followed also 
from the equivariance of the volume form with respect to puUbacks by diffeomor- 
phisms. □ 

5.8. Variation of the covariant derivative. In this section, let V = = V-'^ ^ 
be the Levi-Civita covariant derivative acting on vector fields on M . Since any two 
covariant derivatives on AI differ by a tensor field, the first variation of V-^ ^ is 
tensorial. It is given by the tensor field i'c/,/,)^-^'^ ^ T^iT^M). 

Lemma. The tensor field D{^jj^-^W^ ^ is determined by the following relation hold- 
ing for vector fields X, Y, Z on M : 

Proof. The defining formula for the covariant derivative is 

giVxY, Z) = i [xg{Y, Z) + Yg{Z, X) - Zg{X, Y) 

- g{X, [F, Z]) + giY, [Z, X]) + g{Z, [X, Y]) 
Taking the derivative D(^f f^) yields 

{D^fj,^g){\/xY, Z) + g{{D^fj^)\/){X, Y), Z) 

= 

- {D(fj^)g){X, [Y,Z]) + {D^sMg){Y, [Z,X]) + {D(fj^)g)(Z, [X,Y]) 

Then the result follows by replacing all Lie brackets in the above formula by co- 
variant derivatives using \X, Y] — V xY — VyX and by expanding all terms of the 
form using 

X{{DifM9){y,Z)) = 

{\/xD^fj^^g){Y,Z) + {D(fj^)g){^xY,Z) + {D^jj^)g){Y,V xZ). □ 

5.9. Variation of the Laplacian. The Laplacian as defined in section [3. 11 [ can be 

seen as a smooth section of the bundle L(rimm;nmm) over Imm since for every 
/ G Imm it is a mapping 

A-'^ ^ : Tjlmm — T^Imm. 

The right way to define a first variation is to use the covariant derivative defined 
in section 14.21 



Lemma. For A G r(i(TImm; TImm)), / G Imm and ft, h G T^Imm one hi 
(V, A)(,.) . T,. - V(,.,„,„^^„,,,,.,„„ ,^,.,) 

+ V'(B5(/„r/)/i) -Tr» (fi'l/i.T/jV/i). 
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Proof. Let / be a curve of immersions and h a vector field along /. One has 

A : Imm — > L{Tlmm; Tlmm), Ao/ = A^*» :M^Imm^L(ri mm; TImm). 
Using property |3.7| 5 one gets 



iVfA)(.h) = {Va,iAo f)){h)^ V 9^ Ah- AVg^h 
= -Va,Tr9(v2;i)-AVa,/i 

= Tr {g-\D^fj^^g)g-^V'h) - Ty^iVgy^h) ~ A{Ve^h). 

The term Tr^ (Va^ V^/i) will be treated further. Let X, Y be vector fields on M that 
are constant in time. When they are seen as vector fields on M x M then Vg^X = 
Va^y — 0. Using the formulas from section 3.8 to swap covariant derivatives one 
gets 

{Voy^h){X, Y) = Va, {VxVYh ~ V^^yh) 

= Vx^a^yh + R^{ft,Tf.X)Vyh - Vgyv^yh 
= VxVyVa,h + Vx {R^{ft,Tf.Y)h) + R3{ft,Tf.X)Vyh 
- "^v^YVa.h - "^la^.v^nh - R^{ft,Tf.VxY)h. 

The Lie bracket is 

[dt,VpY] = iD^fj^y){X,Y) 
since (now without the slight abuse of notation) 

^{Ou,{DifMV){X,Y)). 

Therefore 

{We,^^h)(X,Y) = 

= {W^W9MX>y) + ^x{R'{ft,Tf.Y)h)+Ra{ft,Tf.X)Wyh 

- Vp,, ,^,v)(x,Y)/i - R^UuTJ.VxY)h 

= {S/^VqMX. Y) + {VTf.xR')UuTf.Y)h + R^{VxfuTf.Y)h 
+ R^{fuVxTf.Y)h + R^{fuTf.Y)Vxh + R^{fuT f .X)Vyh 

- '^(D,,j^,v)(x.Y)h - R^Uu Tf.VxY)h 

= {W^Wo,h){X, Y) + {WTf.xR'){ft,Tf.Y)h + RS{Wxft,Tf.Y)h 

+ R^{ft, {VTf){X, Y))h + R^{ft,Tf.Y)Vxh + R^UuTf.X)Vyh 

= {V^VdMX.Y) + Vx{R'{fuTf.Y)h)+RS{ft,T}.X)Vyh 

~^(^(/,/t)V)(X,Y)^ 
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Putting together all terms one obtains 

[VfAm - Tr {g-\D^fj^^g)g-'V'h) - Tr^' (v {R^ft,Tf)h) 

^ Tr9 (i?^(/t,r/)V/i) + VTr«(D,,,,^)V)/i 
= Tr {g-\D^fj^^g)g-'V'h) + V* {R^{ft,Tf)h) 

^ Tr9 (R^ifuTf)Vh)+VT,HD,fj^,v)h. 



It remains to calculate Tr®(_D(jj^-) V). Using the variational formula for V from 
section |5.8| one gets for any vector field Z and a (7-orthonormal frame Si 



5(Trni?(/,/.)V),^) 

= ^ ^{VD(^fj^)g) (si (g) Si (E) Z + (g) (g) Z ~ Z (g) 

i 

= -(V*(i?(/,/.)5))(^) - I Tr^(Vzi?y,/,)g) 
= -(V*(^(/,/.)5))(^)-^VzTr«P(/,/,)3) 

= -(v*P(/,/.)5) + 2'^Tr^(^(/J0.9))(^) 
= -.9((V*(i?(/j,)5) + ^dTr^P(/,/,)5))*,^^ 



Therefore 



Tr''(2?(/,/,)V) = -(V*(Z?(/,/,).g) + idTr«(D(/,/,)5))". □ 



6. SOBOLEV-TYPE METRICS 

Assumption. Let P be a smooth section of the bundle L(TImm; Tlmm) over Imm 
such that at every f £ Imm the operator 

Pf : T^Imm — )■ Tjlmm 

is an elliptic pseudo differential operator that is symmetric and positive with respect 
to the H'^ -metric on Imm, 

H'lih,!,)^ [ g{h,k)vo\{g). 

J M 

Note that an elliptic symmetric operator is self-adjoint by [21] 26.2]. Then P 
induces a metric on the set of immersions, namely 

G^{h,k)= / g{Pfh,k)vo\{g) for / G Imm, h,kGTflmm. 
J M 

The metric is positive definite since P is assumed to be positive with respect 
to the TJ^-metric. In this section, the geodesic equation on Imm and Bi for the 
G^-metric will be calculated in terms of the operator P and it will be proven that 
it is well-posed under some assumptions. 
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6.1. Invariance of P under reparametrizations. 

Assumption. It will be assumed that P is invariant under the action of the 
reparametrization group Diff(Af) acting on Imm(M, iV), i.e. 

P = {r^PyP for all G Diff (M). 

For any / G Imm and ip G Diff(M) this means 

Pf^{Tfr^)-'oPf,^oTfr^. 
Applied to h G Tylmni tliis means 

Pf{h) oip = Pfo^iho ip). 

The invariance of P imphes that the induced metric is invariant under the 
action of Diff (M), too. Therefore it induces a unique metric on Bi as explained in 
section 14.71 



6.2. The adjoint of VP. The following construction is needed to express the 
metric gradient H which is part of the geodesic equation. Hf arises from the 
metric G/ by differentiating it with respect to its foot point / G Imm. Since G is 
defined via the operator P, one also needs to differentiate Pf with respect to its 
foot point. As for the metric, this is accomplished by the covariant derivate. For 
P G r(i(TImm; Tlmm)) and m G TImm one has 

V™P G r(i(rimm; TImm)), VP G T{L{THmm; TImm)). 

See section l42l for more details. 

Assumption. It is assumed that there exists a smooth adjoint 

Adj(VP) G r(L2 (TImm; TImm)) 
of VP in the following sense: 

g((V™P)/i,fc)vol(.g) - / g(TO,Adj(VP)(/i,fc))vol(5). 

M JM 

The existence of the adjoint needs to be checked in each specific example, usually 
by partial integration. For the operator P = 1 + AA^, the existence of the adjoint 



will be proven and explicit formulas will be calculated in sections 8.2 and|8.3| 



Lemma. // the adjoint of VP exists, then its tangential part is determined by the 
invariance of P with respect to reparametrizations: 



Adj(VP)(/i, fc)^ = {gi^Ph, k) - g{Wh, Pk)) 
for f G Imm, h,k £ Tyimm 



grad^ giPh, k) ~ {g{PK Vfc) + ^(V/i, Pfc))" 
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= Va,|„ {Pf^pix {h o Fl^)) - Pf^pix (Vs,|„ [h o Flf)) 
= ^a.loiPfih) o Flf) - Pf{Vs,\,{h o Flf)) 
= \/x{Pf{h))~Pf{\/xh) 



Therefore one has for m,h,k € T/Imm that 



/ g{m^,AdiiS/P){h,k)^)vo\ig)^ f 5 (T/.rrJ, Adj (VP) (/i, fc))vol(g) 

= I 9{{VTf.mrP)Kk)Y0\{g) 
JM 

= f g{V^T{Ph) - P{V^rh),k)yol{g) 

= I {g{y,nrPh,k)~giV^rh,Pk))Yo\ig) 
Jm 

= / g(m^,{g{VPh,k)~giVh,Pk))^)Yo\{g). □ 
Jm ^ ' 



6.3. Metric gradients. As explained in section |4.4[ the geodesic equation can 
be expressed in terms of the metric gradients E. and K. These gradients will be 
computed now. We shall use that P/ is invertible on the space of smooth sections. 
This follows because P/ is an elliptic, self-adjoint, and positive operator, see the 
beginning of the proof of theorem |6.6| for a detailed argument. 

Lemma. // Adj(VP) exists, then also H and K exist and are given by 

Kfih, m) = PJ^ [iVmP)h + TrS (g(Vm, T/)) .P/i) 
Hf{h, k) = P7I ( Adj(VP)(/^, k)^ - Tf. {g{Ph, Vfc) + g{Wh, Pk)f 
-^g{Ph,k).TT^{S)). 
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Proof. For vector fields m,h,k on Inim one has 

{y,r.GP){h,k) = Df^f^m) I g{Ph,k)vo\{g) 

J M 

g{P{V„,h),k)vo\{g)- [ g(P/i, V„fc)vol(.g) 



D[f,ni)g{Ph, fc)vol(.g) + / g{Ph, k)D(^j^^)\o\{g) 

M JM 

g{P{Vrah),k)Yo\{g)~ [ g{Ph,V ^k)Yo\{g) 

M JM 

g{V,n{Ph),k)vo\ig)+ [ g{Ph,V mk)yo\{g) 

M JM 



M 



g{Ph, fc)L»(/,„)Vol(g) 



g{P{y„^h),k)Y0\{g)- / g(P/i, V,„fc)vol(5) 

A[ JM 

= [ g((V™P)/i,fc)vol(5)+ / g{Ph,k)D(f^„,)Yo\{g) 

JM JM 



One immediately gets the X-gradient by plugging in the variational formula 5.7 for 
the volume form: 

Kf{h,m) = F^7i((V,„F)/i + Trf {g{Vm,Tf)).Ph 

To calculate the _ff-gradient, one rewrites equation ([T]) using the definition of the 
adjoint: 

(V„G^)(/i,A;)= / g(m,Adj(VP)(^A;))vol(,g)+ / g(P/i, fc)i?(/,„)Vol(<?). 

JM JM 

Now the second summand is treated further using again the variational formula of 
the volume density from section [5 .71 



g(P/i,fc)I?(/.„)V0l(.g) = / g{Ph,k)TT^g{\/m,Tf))vo\{g) 

M JM 

= f g{Ph,k)Tya {Vg{m,Tf) -g{m,VTf))Y0\[g) 

JM 

= ( giPh,k)(^-V*gim,Tf)~g{m,Tr\S))yolig) 



gl {VgiPK k),gim, T/)) vol(5) - / g{PK k)g{m, Tr^(5)) vol(5) 

M JM 

[ g{m, -Tf. grad^giPh, k) - g{Ph, k) Tr^ (5)) vol(g) 

JM 



IM 

Collecting terms one gets that 



G'fiHfih,k),m)^iV^Gnih,k) 

= / g{m,Adi{VP){h,k)-Tf.grad^g{Ph,k)-g{Ph,k)Ti-SiS))vo\{g) 

JM 
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Thus the iJ-gradicnt is given by 

Hf(h, k) = ( M]{VP){h, k) - Tf. grad» g{Ph, k) ~ g{Ph, k). Tr9(5)^ 

The highest order term grad^ g{Ph, k) cancels out when taking into account the 
formula for the tangential part of the adjoint from section [6^ 

Hf{h, k) = ( Adj(VP)(/i, k)^ - Tf. {g{Ph, Vfc) + g{Wh, Pk))^ 

-g{Ph,k).Tr3{S)). □ 



6.4. Geodesic equation on immersions. The geodesic equation for a general 



metric on Imm(M, N) has been calculated in section 4.4 and reads as 



'^oJt = lHf{ftJt)-KfiftJt). 

Plugging in the formulas for H, K derived in the last section yields the following 
theorem. 

Theorem. The geodesic equation for a Sobolev-type metric on immersions is 
given by 

^oJt ^\p-\k&]{VP){fu ft)^ - 2.Tf.g{Pft,Vft)^ ~ g{Pfuft).Tr^iS)) 
- p-'[{Vf,P)ft + Trf {g{VfuTf))Pft 

6.5. Geodesic equation on immersions in terms of the momentum. The 

geodesic equation in terms of the momentum has been calculated in section [4. 5 1 for 
a general metric on immersions. For a Sobolev-type metric , the momentum 
G^ {ft, ■) takes the form 

P = Pft® vol{g) : R T*Imm 

since all other parts of the metric (namely the integral and g) are constant and can 
be neglected. 

Theorem. The geodesic equation written in terms of the momentum for a Sobolev- 
type metric G^ on Imm is given by: 

P = Pft ® vol(.g) 

Va,P = ^ ( Adj (VP) {fuft)^- 2Tf.giPft ,\/ ft)^ - g{Pft , /* ) Tr^ (S)) vol(<?) 

6.6. Well-posedness of the geodesic equation. It will be proven that the ge- 
odesic equation for a Sobolev-type metric G^ on Imm is well-posed under some 
assumptions on P. These assumptions are satisfied for the operator 1 -|- AA^ con- 
sidered in section [sj It will also be shown that (7r,exp) is a diffeomorphism from a 
neighbourhood of the zero section in TImm to a neighbourhood of the diagonal in 
Imm X Imm. 

Before we can state the theorem, we have to introduce Sobolev completions of 
the relevant spaces of mappings. More information can be found in |21j . [9], and in 
[5]. We consider Sobolev completions of r(i?), where — >■ M is a vector bundle. 
First we choose a fixed (background) Riemannian metric g on M and its covariant 
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derivative V^^. We equip E with a (background) fiber Ricmannian metric and 
a compatible covariant derivative V^. Then the Sobolev space H^{E) is the Hilbert 
space completion of the space of smooth sections T{E) in the Sobolev norm 

l|/^llfe = E / (5^®5°)((V^)-''^(V^F/i)vol(5). 

This Sobolev space does not depend on the choices of V*^, and since M 
is compact: The resulting norms are equivalent. 

We shall need the following results (see j8], e.g.): 

Sobolev lemma. If k > dim(Af)/2 then the identy on T{E) extends to a injective 
bounded linear mapping H^'^'p{E) C^{E) where C^iE) carries the supremum 
norm of all derivatives up to order p. 

Module property of Sobolev spaces. If k > dim(Af)/2 then pointwise eval- 
uation II^{L{E,E)) X II^[E) — > II^{E) is hounded bilinear. Likewise all other 
pointwise contraction operations are multilinear hounded operations. 

This allows us to define Sobolev completions of Imm and TImni. In the canonical 
charts for Imm(M, A^) centered at an immersion /o, every immersion corresponds 
to a section of the vector bundle f^TN over M (see [TH section 42]). The smooth 
Sobolev manifold Imm'^(A/, iV) (for k > dim(M)/2 + 1) is constructed by gluing 
together the Sobolev completions H'^{fQTN) of each canonical chart. One has 

Imm'=+^ (A/, N) C Imm'' (Af, N) , f] Imm'' ( Af, N) = Imm( A/, N) . 

k 

Similarly, Sobolev completions of the space Tlmm C C°° {M, TN) are defined as 
i/'^-mappings from M into TN, i.e. Tlmm'' = H''{M,TN). 

Assumption 1. P. WP anrf Adj(VP)^ are smooth sections of the bundles 
L(TImm; Tlmm) (Tlmm; Tlmm) (Tlmm; Tlmm) 

Y Y " 

Imm Imm Imm, 

respectively. Viewed locally in trivializations of these bundles, 

Pfh, {VP)f{h,k), (Adj(VP)/(/i,fc))^ 

are pseudo-differential operators of order 2p in h, k separately. As mappings in the 
footpoint f they are non-linear, and it is assumed that they are a composition of 
operators of the following type: 

(a) Local operators of order I < 2p, i.e., nonlinear differential operators 
A{f){x) = A{x, Vf{x),V'-^f{x), ijf{x),f{x)), 

(h) Linear pseudo- differential operators of degrees U, 
such that the total (top) order of the composition is < 2p. 
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Assumption 2. For each f G Imm(A/, A'^), the operator Pf is an elliptic pseudo- 
differential operator of order 2p for p > which is positive and symmetric with 
respect to the H'^ -metric on Imm, i.e. 

/ g{Pfh,k)Yo\{g)^ I g{h, Pfk)yo\{g) for h, k e Tflmm. 



Assumption 3. P is invariant under the action o/Diff(M). (See section 6.1 for 
the definition of invariance.) 

Theorem. Let p>l and k > dim(Af)/2 + 1, and let P satisfy assumptions 1-3. 



Then the initial value problem for the geodesic equation (6.4) has unique local 
solutions in the Sobolev manifold Imm'^^^^ of H'^'^'^p -immersions. The solutions 
depend smoothly on t and on the initial conditions /(O, . ) and /t(0, . ). The 
domain of existence (in t) is uniform in k and thus this also holds in Imni( A/, N) . 

Moreover, in each Sobolev completion Imm'^^^^, the Riemannian exponential 
mapping exp^ exists and is smooth on a neighborhood of the zero section in the 
tangent bundle, and (7r,exp^) is a diffeomorphism from a (smaller) neigbourhood 
of the zero section to a neighborhood of the diagonal in Imni'^^^'' x Imm'°^^^. All 
these neighborhoods are uniform in k > dim(A/)/2 + 1 and can be chosen H^o+'^p - 
open, for ko > dini(M)/2 + 1. Thus both properties of the exponential mapping 
continue to hold in Iinni(M, N) . 



This proof is partly an adaptation of [201 section 4.3]. It works in three steps: 
First, the geodesic equation is formulated as the flow equation of a smooth vector 
field on each Sobolev completion Tlmm*^^^^. Thus one gets local existence and 
uniqueness of solutions. Second, it is shown that the time-interval where a solution 
exists does not depend on the order of the Sobolev space of immersions. Thus one 
gets solutions on the intersection of all Sobolev spaces, which is the space of smooth 
immersions. Third, a general argument involving the inverse function theorem on 
Banach spaces proves the claims about the exponential map. 



Proof. By assumption 1 the mapping Pfh is of order 2p in / and in h where / is 
the footpoint of h. Therefore / i— > extends to a smooth section of the smooth 
Sobolev bundle 

L(rimm*=+2p.yjj^^fc I inini'^+2p) ^ Imm'^+^P , 

where Tlmm'^ | Imm*''*''^^ denotes the space of all H'' tangent vectors with foot 
point a H^+'^P immersion, i.e., the restriction of the bundle Tlmm'^ — >■ Imm'^ to 
Imm''+^'' C Imm''^ 

This means that Pf is a bounded linear operator 

Pf e L{H''+^P{f*TN),H''{f*TN)) for / e Imm'^+^f. 

It is injective since it is positive. As an elliptic operator, it is an unbounded operator 
on the Hilbert completion of T/Imm with respect to the iJ'^-metric, and a Fredholm 
operator H^+'^p — > for each k. It is selfadjoint elliptic, thus by [HI theorem 
26.2] it has vanishing index. Since it is injective, it is thus also surjective. 
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By the implicit function theorem on Banach spaces, / i— > is then a smooth 
section of the smooth Sobolev bundle 

L(TImm'^ | Imm'^^^^; Tlmm'^^^^) — )■ Imm''^^'' 

As an inverse of an elliptic pseudodifFerential operator, Py^ is also an elliptic 
pseudo-differential operator of order —2p. 

By assumption 1 again, {\/P)f{m, h) and ( Adj(VP)/(m, h))^ are of order 2p in 
/, m, h (locally). Therefore f ^ Pf and / i-> Adj(VP)-'~ extend to smooth sections 
of the smooth Sobolev bundle 

L^{TlmYn^+^P-Tlmm'' \ Imm'^'+^P) ^ Imm'^+^f 

Using the module property of Sobolev spaces and counting the order of all re- 
maining terms in the geodesic equation |6.4[ one obtains that the Christoffel symbols 

]^Hf{h, h) - Kf{h, h) = ^P-I ( Adj(VP)(/i, h)^ - 2.Tf.g{Ph, Vhf 

- g{Ph, h).Tr3{S) ~ {WhP)h - Tr^ {g{Wh, Tf ))Ph) 

extend to a smooth (C°°) section of the smooth Sobolev bundle 
L2y^^(rimm'=+2P;TImm^+2P) ^ Imm'^+^f 

Thus h ^ \H{h,h) - K{h,h) is a smooth quadratic mapping TImm TImm 
which extends to smooth quadratic mappings Tlmm*^^^^ — ^ Tlmm'^"'"^^ for each 
k > '^'"''^^ + 1. The geodesic equation 

^lft^lHf{fuft)-Kf{ftJt) 

can be reformulated using the linear connection : TA^ x jv TN — >■ TTN (hori- 
zontal lift mapping) of V^, see [T7J section 24.2]: 

dtft = c(^^Hf{ftJt)-Kf{ftJt),ft). 

The right-hand side is a smooth vector field on Tlmni'"'^^^, the geodesic spray. Note 
that the restriction to Tlmm'"'^^^^^ of the geodesic spray on Tlmm*^^^^ equals the 
geodesic spray there. By the theory of smooth ODE's on Banach spaces, the flow 
of this vector field exists in Plmm'^^^^ and is smooth in time and in the initial 
condition. 

Consider a C°° initial condition Hq G TImm with foot point /o G Imm. Sup- 
pose the trajectory F\^{ho) of geodesic spray through these initial conditions in 
Tlmm'^'*'^'' maximally exists for t € (—0^,6^), and the trajectory Fl(^^(ft,o) hi 
Tlmm*^^^^^^ maximally exists for t G {—ak+i,bk+i) with b^+i < bk, say. By 
uniqueness of solutions one has Fl(^^(/io) = Fl^(/io) for t G (— flfc+i.f'fc+i). We now 
write V for the covariant derivative induced by g on A^ and the background metric 
g on M. Let X be a vector field on M. Applying —: V x to the geodesic 
equation and swapping covariant derivatives yields: 

(A) - V|(ii//(,A,/0 - Kfifuft)) + RHfuTf.X)ft 
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Note that ixVo,VS/t = Vg^V^/* - V^^ ^ft = V|^V^/( - 0. Thus we can omit 

X and rewrite (A) as an equation for Tf. We aim to rewrite equation (A) as a 
linear first order equation for the highest derivative whose coefficients arc given 
by Flj(/io) and thus exist beyond (— afe+i, 6fe+i). For this we have to pass to one 
(ore more) canonical chart for Imm and the induced trivializations of all bundles 
as before and in assumption (1). Then / itself has values in a vector space and we 
may regard (A) as a vector valued 1-form on M. So we rewrite (A) as: 

Va,VVt = V'{\Hj{fuft) - Kj(Jufi)) + R'{.ft,Tf),ft 

(B) = mv'^Tf) + Y^if, ftKV^P+'ft) + Y^f, h). 
We claim that (B) consists of: 

• The smooth expression (/, ft)iy^^^'^ f) which is linear and of order in V^^"*"^/ 
and where Y^{f, ft) is of order < 2p in /, ft] order here means that the expression 
prolongs continuously to the corresponding Sobolev spaces. 

• The smooth expression Y^{f, ft){\7^P^^ ft) which is linear and of order in 
V^P+Vt and where Y'^{f, ft) is of order < 2p in /, ft. 

• The smooth expression Y^{f, ft) of order < 2p in /, ft- 

To see this we claim that the highest derivatives of order 2p + 1 of / and ft appear 
only linearly in (A). This claim follows from assumption i: 

(a) For a local operator we can apply the chain rule: The highest derivative of / 
appears only linearly. 

(b) For a linear pseudo differential operator A of order k the commutator [V, A] is 
a pseudo-differential operator of order k again. 

On the left hand side of (B) we write Tf =: u and V/f =: v. On the right hand 
side of (B) we write V^^T/ = V^p^ and V^^'+Vt = V^^'t; for the highest derivatives 
only. Then the system (B) becomes: 

Vg^U = V 

(C) Ve,v = Y\f,ft){V'^u) + Y^f,ft){V'Pv) + Y^f,ft). 

The coefficients /, ft in (C) exist for t G (— a^, bk) as Flf (/iq). Then (C) is a bounded 
and smooth inhomogeneous linear ODE for {u,v) G n^iM,Tlmm''+^P), I.e., m a 
Banach space. This equation therefore has a solution {u{t), v{t)) for all t for which 
the coefficients exists, thus for all t G (a^, bk), which is unique for the initial values 
Up = T/o and v„ = Who- The limit Mmt/'b.+i {u{t),v{t)) exists in n^M, Tlmm'^+^P) 
and by continuity it equals (VTf, 'V ft) for t — bk+i- Thus the flow line Flj (/iq) was 
not maximal and c;an be continued. So assuming bk+i < bk leads to a contradiction, 
and thus {—ak+i,bk+i) = (—0^,6^.). Iterating this procedure one concludes that 
the flow line FY^{ha) exists in nfc> d.„(M) , ^ rimm'=+^^' = Tlmm. 

It remains to check the properties of the Riemannian exponential mapping exp^. 
It is given by exp^(/i) = c(l) where c{t) is the geodesic emanating from value / with 

initial velocity h. Let ko > dim(M)/2 + 1 and k > ko. On each space Tlmm*^"*"^^, 
the properties claimed follow from local existence and uniqueness of solutions to 
the flow equation of the geodesic spray, from the form of the geodesic equation 
ftt = ^H{ft, ft) — K{ft, ft) when it is written down in a chart, namely linearity in 
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ftt and bilinearity in ft , and from the inverse function theorem which holds on each 
of the Sobolev spaces Imm'^"''^^. See for example [T71 22.6 and 22.7,] for a detailed 
proof which works without any change in notation. 

Imm''°+^P contains Imm'^^+^P for k > ko. Since the spray on Imm'^""''^'' restricts 
to the spray on each Imm'"^^^, the exponential mapping exp^ and the inverse 
(tt, exp-'^)"^ on Imm'^^^^^ restrict to the corresponding mappings on each Imm''^^''. 
Thus the neighborhoods of existence are uniform in k and can be chosen H^o+'^P- 
open. □ 



6.7. Momentum mappings. Recall that by assumption, the operator P is in- 
variant under the action of the reparametrization group Diff(Af). Therefore the 
induced metric is invariant under this group action, too. According to [201 
section 2.5] one gets: 

Theorem. The reparametrization momentum, which is the momentum mapping 
corresponding to the action o/ Diff(Af) on Imm(M, iV), is conserved along any 
geodesic f in Imm(M, iV); 

yX e X(M) : / 5(T/.X,P/0vol(5) 

or equivalently 

g{iPftV)vo\{g) e T{T*M vo1(M)) 

is constant along f. 



6.8. Horizontal bundle. The splitting of TImm into horizontal and verti cal sub- 
spaces will be calculated for Sobolev-type metrics . See section 4.7 for the 
definitions. By definition, a tangent vector h to f E Imm(M, N) is horizontal if 
and only if it is G^-perpendicular to the Diff (M)-orbits. This is the case if and 
only if g{Pfh{x),Txf.Xx) = at every point x S M. Therefore the horizontal 
bundle at the point / equals 

{h e T/Imm : Pfh{x)±T^f{T^M) for ah x e M} ^ {h : {PfhY = O}. 

Note that the horizontal bundle consists of vector fields that are normal to / when 
P — Id, i.e. for the iJ^-metric on Imm. 

Let us work out the G^-decomposition of h into vertical and horizontal parts. 
This decomposition is written as 

h = Tf.h""' + h^°\ 

Then 

Pfh = PfiTf.h''") + Pfh^'°' with {Pfh^ = {Pj{Tf.h^")Y + 0. 
Thus one considers the operators 

Pj : X{M) ^ X(Af), Pj{X) = {PfiTf.X)f, 

Pf,i_ : X(Af) ^ r(Nor(/)) C G^{M,TN), PfAX) = {PiiTf.X))^. 
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The operator is unbounded, positive and symmetric on the Hilbert completion 
of T^Imm with respect to the _ff'^-metric since one has 

/ g(PjX,Y)vo\{g)^ f g{Tf.PjX,Tf.Y)vo\{g) 
= f g{Pf{Tf.X),Tf.Y)yol{g) 
= / g{PjY,X)vol{g), 

JM 

[ g{PjX,X)vo\{g)^ [ g{Pf{Tf.X),Tf.X)vo\{g)>0 ifX^O. 
Jm Jm 

Let a^f and a^f denote the principal symbols of Pf and Pj , respectively. Take 
any x & M and ^ G T*M \ {0}. Then cr^^ {£,) is symmetric, positive definite on 
{Tf^x)N,g). This means that one has for any h,k (z Tf^^-fN that 

g{aPf{i)h,k)^g{h,aPf{Ok), g{a'''{Oh,h) > for /i ^ 0. 

The principal symbols a^f and a^f are related by 

g{aPl iOX.Y) ^giTf.a^I {OX,Tf.Y) ^g{<7Pf{0Tf.X,Tf.Y), 

where X,Y ^ T^M. Thus a^i (^) is symmetric, positive definite on [T^M^g). 
Therefore Pj is again elliptic, thus it is selfadjoint, so its index (as operator 
jjk+2p _^ ^fe-j vanishes. It is injective (since positive) with vanishing index (since 
self-adjoint elliptic, by |21] theorem 26.2]) hence it is bijective and thus invertible 
by the open mapping theorem. Thus it has been proven: 

Lemma. The decomposition oj h ^ TImm into its vertical and horizontal compo- 
nents is given by 

h^^^^{pj)-\{Pfhy), 

f^hor ^ ^ _ Tf.h""' =h- Tf.{Pj)-^{[PfhY) . 



6.9. Horizontal curves. To establish the one-to-one correspondence between hor- 
izontal curves in Imm and curves in shape space that has been described in theo- 
rem 



4.8 one needs the following property: 



Lemma. For any smooth path f in Imm(M, N) there exists a smooth path if in 
Diff(M) with ip{0, . ) — Um depending smoothly on f such that the path f given 
by f{t,x) — f{t,(p(t,x)) is horizontal: 

g{Pf{dtf),Tf.TM)=0. 

Thus any path in shape space can be lifted to a horizontal path of immersions. 



The basic idea is to write the path ip as the integral curve of a time dependent 
vector field. This method is called the Moser-trick (see [TBI Section 2.5]). 

Proof. Since P is invariant, one has {r'^)*P = P or Pfo^{u o ip) — (Pfu) o ip for 
(fi S Diff(A/). In the following foip will denote the map f{t, (p{t, x)), etc. One looks 



34 



MARTIN BAUER, PHILIPP HARMS, PETER W. MICHOR 



for ip as the integral curve of a time dependent vector field S,{t^x) on M, given by 
'^t — S, ° 'f- The following expression must vanish for all a; e M and Xj. e T^M: 

= g[Pfo^{dt{f o ^)){x),T{f o (^).X,) 

^g[Pfo^{{dtf) o ^ + Tf.{dtv)){x),T{f o ^).X.,) 

= 9{{{Pfidtf)) + PfiTf.O) {v{x)),Tf o T^.X^) 

Since Ttp is surjective, Tip.X exhausts the tangent space T^i^^-^M, and one has 

{{Pf{dJ))+Pf{Tf.O){^{x)) ± /. 

This holds for all x E M, and by the surjectivity of (p, one also has that 

{(Pf{dJ)) + Pf{Tf.O){x) ± f 

at all X G M. This means that the tangential part (Pf{dtf) + Pf{Tf.^))^ vanishes. 
Using the time dependent vector field 

and its flow ip achieves this. □ 

6.f0. Geodesic equation on shape space. By the previous section and theo- 
rem |4.8[ geodesies in Bi correspond exactly to horizontal geodesies in Imm. The 
equations for horizontal geodesies in the space of immersions have been written 
down in section [4. 9| Here they are specialized to Sobolev-type metrics: 

Theorem. The geodesic equation on shape space for a Sobolev-type metric is 
equivalent to the set of equations 

ft = f}°' e Hor, 
(VaJt)'^°'- - lp-'[Adi{VP){ftJt)^~g{PftJt).Tr'{S)) 

- P''{{i^f.P)ft)^ - Tr^ {g{yfuTf))Pft), 
where f is a horizontal path of immersions. 

These equations are not handable very well since taking the horizontal part 
of a vector to Imm involves inverting an elliptic pseudo-differential operator, see 
section 16.81 However, the formulation in the next section is much better. 



6.11. Geodesic equation on shape space in terms of the momentum. The 



geodesic equation in terms of the momentum has been derived in section 4.10 for a 
general metric on shape space. Now it is specialized to Sobolev-type metrics using 
the formula for the 77-gradient from section |6.3[ 



As in section 6.5 the momentum G {ft, •) is identified with Pft vo\{g). By 
definition, the momentum is horizontal if it annihilates all vertical vectors. This is 
the case if and only if Pft is normal to /. Thus the splitting of the momentum in 
horizontal and vertical parts is given by 

Pft (g> vol(<?) = (Pft)^ vol(g) + Tf.[PftV vol(<?). 
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This is much simpler than the sphtting of the velocity in horizontal and vertical 
parts where a pseudo-differential operator has to be inverted, see section [678] Thus 
the following version of the geodesic equation on shape space is the easiest to solve. 

Theorem. The geodesic equation on shape space is equivalent to the set of equations 
for a path of immersions f : 

p = Pft®Yo\{g), Pft^{Pf,)^, 
(Va.p)'^°^ = l^AdjiV P){ftJt)^-giP ft Jt).Ti-'{S)) ®vol(g). 

The equation for geodesies on Imm without the horizontality condition is 

Va,p = ^ ( Adj(VP)(/,, /i)^ - 2Tf.giPft,Vft f - g{PftJt).Tv^iS)) ® volCg), 

see section |6.5| It has been proven in section |4.9| that the vertical part of this 
equation is satisfied automatically when the geodesic is horizontal. Nevertheless 
this will be checked by hand because the proof is much simpler here than in the 
general case. 

If ft is horizontal then by definition Pft is normal to /. Thus one has for any 
X e X{M) that 

g{{Vo,PftV,X)=g{Wa,Pft,Tf.X) = 0-g{Pft,Wa,Tf.X) 

= ^9iPft,'^xft) = ^g{9iPft,Vft)Kx). 

Thus 

(Va,p)™'* = ((Va,P/t) ® vol(g) + Pft ® D^fj^^yo\{g)y'''' 

= Tf.{Vs,PftV ® vol(.g) + Tf.{Pft)^ (g> D^fj^)Y0\{g) 
= -Tf.g{Pft,Vftf®yo\{g)+Q, 
which is exactly the vertical part of the geodesic equation. 

7. Geodesic distance on shape space 

It came as a big surprise when it was discovered in [19J that the Sobolev metric of 
order zero induces vanishing geodesic distance on shape space Bi. It will be shown 
that this problem can be overcome by using higher order Sobolev metrics. The 
proof of this result is based on bounding the G^-length of a path from below by 
its area swept out. The main result is in section f7.6| The same ideas are contained 
in [U section 2.4], [21 section 7] and [HI section 3]. 

7.1. Geodesic distance on shape space. Geodesic distance on Bi is given by 

dist^'p(^^o,i^i) -infiGp(^)' 

F 

where the infimum is taken over all F : [0, 1] Bi with F(0) = i^o and F{\) = Fi. 
L^p is the length of paths in Bi given by 

L^AF)^J^ ^GP{Ft,Ft)dt ioY F -.[0,1]^ B,. 
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Letting tt : Imm — >■ Bi denote the projection, one has 

r-l 



when / : [0, 1] — ^ Imm is horizontal. In the following sections, conditions on the 
metric ensuring that dist^p separates points in Bi will be developed. 

7.2. Vanishing geodesic distance. 

Theorem. The distance dist^o induced by the Sobolev metric of order zero 
vanishes. Indeed it is possible to connect any two distinct shapes by a path of 
arbitrarily short length. 

This result was first established by Michor and Mumford for the case of planar 
curves in |19j . A more general version can be found in [18j . where the same result 
is proven also on difleomorphism groups. 

7.3. Area swept out. For a path of immersions / seen as a mapping / : [0, 1] x 
M — >• TV one has 

(area swept out by /)= / vol(/(-, = / / \\ft\\^o\(g)dt. 

J[0,l]xAf Jo Jm 

7.4. Area swept out bound. 

Lemma. Let G^ be a Sobolev type metric that is at least as strong as the H^- 
metric, i.e. there is a constant Ci > such that 

||^|IgJ= — ^1 ll'*ll_ffo = I ^^^^^^^'^^^^^ ^ ^ TImm. 

Then one has the area swept out bound for any path of immersions f: 



Ci {area swept out by f) < max Y^Vol(/(i)) 

The proof is an adaptation of the one given in section 7.3] for almost local 
metrics. 



Proof. 



Lp^if)^ f\\ft\\G^dt>C^ fwftWHodt 
Jo Jo 

/Vt^LoC^i-Ci r ( f \\ft^\o\{g)y dt 
Jo Jo ^Jm ' 

>C, [' ( f vol(g))"' / l.\\f,^vo\{g)dt 
'0 ^Jm ' Jm 



> 



Ciniinf / vol(g)) ' / vol(/(., ^s) 

* ^Jm ' J[0,l]xM 

C\{ max / No\{g)\ ^ (area swept out by /) □ 
^ * Jm ' 
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7.5. Lipschitz continuity of vvol. 

Lemma. Let he a Soholev type metric that is at least as strong as the - 
metric, i.e. there is a constant C2 > such that 



\\h\\aP>C2\\h\\^i=C2^J g{{l + A)h,h)vol{g) for all h e Tlmm. 
Then the mapping 

VVoi: (B,,dist^'p) ^M>o 
is Lipschitz continuous, i.e. for all Fq and Fi in Bi one has: 

VVol(Ti) - VVol(Fo) < -^dist^'p(Fo,Fi). 
For the case of planar curves, this has been proven in [501 section 4.7]. 
Proof. 

dtYol - (div3(/7) - Tr^'(5)))vol(5) 

= 0+ / 5(/t,V*T/)vol(.g)= / {9>9){^ft,Tf)vo\{g) 

JM JM 



< 



M V JM 



l|V/t||'o^^vol(ff)W / ||r/||^o^^vol(g) 



< IMhi Vyol<^\\ft\\aP VVol. 



Thus 



By integration one gets 



yvoi(A) - ^voic^y - dtVyoiiDdt < £ ^ ii/*iIg- = ^ iG""(/)- 

Now the infimum over ah paths / : [0, 1] Imni with 7r(/(0)) — Fq and 7r(/(l)) = 
Fi is taken. □ 

7.6. Non-vanishing geodesic distance. Using the estimates proven above and 
the fact that the area swept out separates points at least on B^, one gets the 
following result: 

Theorem. The Soholev type metric induces non-vanishing geodesic distance 
on Be if it is stronger or as strong as the -metric, i.e. if there is a constant 
G > such that 



Hgp > Cll^llffi = cJ J g{{l-\- A)h,h)vo\{g) for all h e Tlmm. 
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Proof. By lemma [7!4| we have 



Ci (area swept out by /) < max y^Vol(/(t)) .L^'°(/). 



Now we use the Lipschitz continuity 7.5 of v Vol and that area swept out separates 



points on Bf.. □ 

8. SOBOLEV METRICS INDUCED BY THE LAPLACE OPERATOR 

The results on non- vanishing geodesic distance from the previous section lead us 
to consider operators P that are induced by the Laplacian operator: 

P =1 + AAP, P £ r(i(TImm; Tlmm)) 

for a constant A > 0. (See section 3.11| for the definition of the Laplacian that is 
used in this work.) At every / e Imm, Pf is a positive, selfadjoint and bijective 
operator of order 2p acting on T/Imm T{f*TN). Note that A depends smoothly 
on the immersion / via the puUback-metric f*g, so that the same is true of P. P is 
invariant under the action of the reparametrization group Diff (A/). It induces the 
Sobolev metric 

J{h,k)^ f g{Pf{h),k)Yolig)^ f g({l + AiAf'^r)h,k)vol{rg). 



'm Jm 
When A = 1 we write HP := 0^+'^" . 

In this section we will calculate explicitly for P = 1 + AA^ the geodesic equation 
and conserved momenta that have been deduced in section [6] for a general operator 
P. The hardest part will be the partial integration needed for the adjoint of VP. As 
a result we will get explicit formulas that are ready to be implemented numerically 



8.1. Other choices for P. Other choices for P are the operator P = 1+j4(V*)^V^ 
corresponding to the metric 

Gf{h, k)^ f {g{h, k) + Ag{VPh, VPk))vo\{g), 

JM 

and other operators that differ only in lower order terms. Since these operators all 
have the same principal symbol, they induce equivalent metrics on each tangent 
space T/Imm. It would be interesting to know if the induced geodesic distances on 
Bi are equivalent as well. 



8.2. Adjoint of VP. To find a formula for the geodesic equation one has to cal- 
culate the adjoint of VP, see section |6.4[ The following calculations at the same 



time show the existence of the adjoint. It has been shown in section 6.2 that the 
invariance of the operator P with respect to reparametrizations determines the 
tangential part of the adjoint: 

Adj(VP)(/i, k)Y = grads g{Ph, k) - {g{Ph, Vfc) + g{Vh, Pk)f . 
It remains to calculate its normal part using the variational formulas from sectionjS] 

In the following calculations there will be terms of the form Tr(g~^sig~^S2), 
where si, S2 are two-forms on M. When the two-forms are seen as mappings TM — >■ 
T*M, they can be composed with g^^ : T*M -> TM. Thus the expression under 
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the trace is a mapping TM — > TM to which the trace can be apphed. When 
one of the two-forms is vector valued, the same tensor components as before are 
contracted. For example when h G r{f*TN) then S2 = W^h is a two-form on M 
with values in f*TN. Then in the expression Tr^g^^ .si.g^^ .S2) only TM and T*M 
components are contracted, whereas the f*TN component remains unaffected. 



g(m^,Adj(VP)(/i,fc))vol(g) = / g{{V^^P)h,k)Yo\ig) 

M JM 
p-1 



^AY^I g((V„,.A)AP-'-i/i,A'fc)vol(g) 

/ 5(Tr(5-^i^(/,™.)5.5-'v2AP-'-i;i),A'fc)vol(g) 

p-1 „ 

/ 5(v*i?^(m^,r/)Af-»-i/i,A*fc)vol(5) 

p-1 „ 

^E / 5(TrS(^^(m-L,T/)VAP-'-i/i),A'fc)vol(g) 
= ^E / Tr(.9-i.i?(^.„.)g..g-ig(V2Af-^-i/i,A*fc))vol(g) 

p-1 „ 

^E / (.9?®ff)(vA?'-^-i/i,(V*%,„.).g)® A'fc)vol(g) 
^E / (5?®5)(vAf-'-i/j,-dTr^p(/,„,.)g)®A^fc)vol(5) 
^E / (.9?®5)(i?>i^,r/)Af-*-i/i,VA'fc)vol(g) 
/ ^(TrS (i?9(m-L,T/)VAP-'-i;i),A^fc)vol(g) 



1=0 

p-1 



1=0 ' 

Using the following symmetry property of the curvature tensor (see |T71 24.4.4]): 
g{K\X, Y)Z, U) = -g{W{Y, X)Z, U) = -g{W{Z, U)Y, X) 

yields: 

' g{m^,Adi{VP){h,k))vo\{g) = 

M 

= ^E//°(^(/.™-).9,5(V'Af-^-i/i,A'A;))vol(.g) 
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p-1 

I M 



-^E / 5?(5(VAf-^-i/i,A^fc),V*D(^,„.)5)vol(5) 

-^E / ff?(5(VAf— i/i,A^fc),JvTr«(i?(/,„.)5))vol(5) 
i=o-^^ ^ ^ ^ 

+ ^E / ^(Trs (i?9(Af-*-i/i,VA*A;)r/),m-L)vol(5) 

-^E / 5(TrS(i?^(VAf-^-i/i,A*A;)r/),m^)vol(5) 

^E / !?2(^(/,m-)5,5(V'A^'— ift,A^fc))vol(5) 
i=o ^ ' 

-^E / 52(v5(VA^'-'-i/.,A^fc),Z)(/,„.)5)vol(ff) 
i=o -^^ ^ ^ 

- ^ E / (v*5(VA^'-^-i/i, A^fc)) Tr«(i?(/,„.)5)vol(5) 
+ AV / ^(TrS (i?9(AP-'-i/i,VA'A;)T/),m-L)vol(5) 
-^E / sfTrS (i?s(VAP-'-i/i,A^A;)T/),m-L)vol(5) 

^E / ff2°(^(/,n.-)5,5(V2Af-'-i/i,A'fc))vol(5) 

-^E / 52°(ff(V'AP-*-ift,A*fc),Z)(/,„.)5)vol(5) 
p-i . 

-^E / 52(fl(VAP-'-i/i,VA^fc),D(/,„.)ff)vol(5) 

- ^ E / (v*5(VAf-^-i/i, A^fc)) Tr^(i^(/,„.)<7)vol(5) 

+ ^E / 5(Tr«(i?^(A^'-'-^/i,VA'fc)T/),m^)vol(5) 
p-i „ 

-^E / Tr-'' (^^(VAP-'-i/i,A'fc)T/),m-L)vol(5) 
-A^ / 5°(%,™-)5,5(VAf-'-i/i,VA'fc))vo%) 
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M 



- 4 E / (v*5(VAf-^-i/i, A'k)) Tr«(D(/,„.)5)vol(ff) 
i=o •'^ 

+ / ^(Trs (i?9(Af-*-i/i,VA'fc)r/),m-L)vol(5) 

p-i . 

-^V / ff(TV«(i?^(VAf-'-i/i,A*fc)r/),m^)vol(5) 
-^E / 92(-mm^,S),g{VAP-'-^h,VA'k))Y0l{g) 

- 2 E (v*5(VAf-*-ift, A^fc)) Tr« ( - 2.5(m^, 5))vol(5) 

p-i „ 

+ ^E / sfTrS (i?9(Af-^-i/i,VA'A:)r/),m-L)vol(5) 

p-i „ 

-AV / ^(TrS (ES(VAf-'-i/i,A*A:)T/),m-^)vol(5) 
/ 5(m^,2AVTr(5-i55-iff(VAf-^-i/i,VA*fc))) 

+ / 5(m^,^E(^*^(^^''"'"'^'^'^))Tr'(^))^oKff) 
■^^ ^ i=o ^ 

+ ^E / (i?9(AP-'-^/i,VA'A;)T/),m-L)vol(5) 

- ^ E / 5 f (-R^(VAf-'-i/i, A*A:)T/) , m-^) vol(5). 



Prom this, one can read off the normal part of the adjoint. Thus one gets: 
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Lemma. The adjoint of VP defined in section \67^ for the operator P = 1 + AA^ 
is 

p-i 

Adj(VP)(/i,fc) = 2AY,TT{9-'Sg-'g{VAP-'-'h,VA'k)) 

i=0 

p-1 



{V*g{VAP-'-^h,A'k))TrS{S) 

1=0 

p-i 

A ^ Tr^ {R^{AP~'~^h, VA'k)Tf) 



i=0 

p-1 



A^Trf {R3(yAP^'^^h,A'k)Tf) 

i=0 

Tf. [ grad'' g{Ph, k) ~ {g{Ph, Vfc) + glVh, Pk))^ 



8.3. Geodesic equations and conserved momentum. The shortest and most 
convenient formulation of the geodesic equation is in terms of the momentum p = 



(1 + AAP)ft ® vol(g), see sections and 6.11 



Theorem. The geodesic equation on Imm(M, A^) for the -metric with P 
1 + AAP is given by: 

p-(l + ^AP)/t®vol(5), 

/ p-i 

Vg,p = U^Tr (5-i5g-ig(V(Af— Vt), VA7t)) 
\ i=0 

A 



+ 9 E (V*5(V(A^-^-Vt), A7,))- Tr^(^) 



i=0 
p-1 



+ 2^^Trf {R^AP-'-'ft,VA^ft)Tf) 



-g{Pfuft) Tr^'(^) - Tf.g{Pfu V ftf ® vol(g). 



This equation is well-posed by theorem \6.6\ since all conditions are satisfied. For 
the special case of plane curves, this agrees with the geodesic equation calculated in 
[201 section 4.2]. 

P — 1 + AAP and consequently are invariant under the action of the 
reparametrization group Diff(Af). According to section 6.7 one gets: 

Theorem. The momentum mapping for the action of Diff (M) on Imm(M, A^) 

^(((1 + AAP)ftf^ ® vol(.g) e T{T*M ®m vo1(M)) 
is constant along any geodesic f in Imm(Af, A^). 



The horizontal geodesic equation for a general metric on Imm has been derived 
in section |4.10[ In section |6.11| it has been shown that this equation takes a very 
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simple form. Now it is possible to write down this equation specifically for the 
operator P 1 + AA^: 

Theorem. The geodesic equation on shape space for the Sobolev-metric with 
P = 1 + AAP is equivalent to the set of equations 



p^Pft®vo\{g), Pf,^(Pf,)^, 



p-1 



(Va.p) 



hor 



A^Tr(5-i55-^5(VAf-^-Vt,VA7,)) 



1=0 



A 



p-1 



^ (V*g(VAf-'- Vt, A7t)) Trf (^) 



2A^Tr5 (i?nA?'-^-V„ VA7t)T/) 



-giPftJt).TrS{S) ®vol(g), 



where f is a path of immersions. For the special case of plane curves, this agrees 
with the geodesic equation calculated in [20', section 4.6] . 



9. Surfaces in n-sPACE 

This section is about the special case where the ambient space N is M". The flat- 
ness of leads to a simplification of the geodesic equation, and the Euclidean mo- 
tion group acting on M" induces additional conserved quantities. The vector space 
structure of M" allows to define a Frechet metric. This metric will be compared to 



Sobolev metrics. Finally in section 9.5 the space of concentric hypcr-sphercs in . 
is briefly investigated. 



9.1. Geodesic equation. The covariant derivative on M" is but the usual 
derivative. Therefore the covariant derivatives Vg^/t ^dtP in the geodesic 
equation can be replaced by fu and pt, respectively. (Note that Imm(M, M") is an 
open subset of the Frechet vector space C°°(M, E").) Also, the curvature terms 
disappear because M" is flat. Any of the formulations of the geodesic equation 
presented so far can thus be adapted to the case N = M" . 

We want to show how the geodesic equation simplifles further under the addi- 
tional assumptions that dim(M) = dim(A^) — 1 and that M is orientable. Then it is 
possible define a unit vector field to M. The condition that ft is horizontal then 
simplifies to Pft = a.v for a € C°°{M). The geodesic equation can then be written 
as an equation for a. However, the equation is slightly simpler when it is written 
as an equation for a.vo\{g). In practise, vol(g) can be treated as a function on M 
because one can identify vol(5) with its density with respect to du^ A ... A du"'~^, 
where (u^, . . . , it"~^) is a chart on M. Thus multiplication by vol(g) does not pose 
a problem. 
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Theorem. The geodesic equation for a Sobolev-type metric on shape space 
_Bi(Af, ]R") with dim(Af) = n — 1 is equivalent to the set of equations 

Pft = a.iy 

dt{a.vo\{g)) = { AdiiV P)i ft, ft), ^) - \giPft,ft)9{Tr'iS),iy), 
where f is a path in Imni(M, M") and a is a time- dependent function on M. 



Proof. Applying g{-,v) to the geodesic equation 6.11 on shape space in terms of 
the momentum one gets 

dt{a.vo\ig)) = dt g{Pft vol(.g), i^) 

= 5(Va, {Pft ® vol(5)), + g{Pft (g> vol(5), Va.z^) 

= ^g{Ad]iVP){fuft),:^) - lg{Pft,ft)g{TT^iS),:y)+0. □ 

Let us spell this equation out in even more details for the iJ^-metric. This is the 
case of interest for the numerical examples in section |11[ 

Theorem. The geodesic equation on shape space i?i(A/, M") for the Sobolev-metric 
G^ with P = 1 + AA is equivalent to the set of equations 

Pft = a.v 

dt{a.^o\{g)) = (Agl{s,g{V fuV ft)) - ^ ( ll/t ll| + ^ || V/* ||^„^^ )) vol(5), 

where f is a path of immersions, a is a time- dependent function on M , s = g{S, v) £ 
T{T2M) is the shape operator, L — g~^s G T{TiM) is the Weingarten mapping, 
and Tr(L) is the mean curvature. 



Proof. The fastest way to get to this equation is to apply g{-,v) to the geodesic 



equation on Imm from section 8.3 This yields 

dt{a.Mo\{g)) = (ATr(.g-i.s..g-ig(V/t,V/t)) + | (V*5(V/t, /*)) Tr(L) 
-\g{Pft,ft).ML)yo\{g) 
Ag°{s,g{Vft,Vft)) I Tr^ (ff(VVt,/t)) Tr(L) 
- ^Tr^ (5(V/t,V/0)Tr(i)- ^5((l + AA)/,,/0Tr(i))vo%) 
= (AgO(s,g(V/t, V/0) - I Tr^ (g(V/t, V/*)) Tr(L) 
-^5(/t,/t)-Tr(i))vol(g). 
Notice that the second order derivatives of ft have canceled out. □ 
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9.2. Additional conserved momenta. If P is invariant under tlie action of tlie 
Euclidean motion group M" xi SO(n), then also the metric is in invariant under 
this group action and one gets additional conserved quantities as described in [201 
section 2.5]: 

Theorem. For an operator P that is invariant under the action of the Euclidean 
motion group M" xi SO(n), the linear momentum 



and the angular momentum 



f PftV0l{g)e{R"r 
Jm 

VXeso(n) : / 5(X/,P/t)vol(5) 
Jm 



or equivalently / (/A P/f)vol(.g) G /\ M" = so{n)* 

JM 

are constant along any geodesic f in Imm(Af, M"). The operator P — 1 + AA^ 
satisfies this property. 

9.3. Prechet distance and Finsler metric. The Frechet distance on shape space 
B,(M,M") is defined as 

distf;(Fo,Fi)- inf Ifo-fih^, 

Joji 

where the infimum is taken over all /o, fi with 7r(/o) = Fq, 7r(/i) = Fi. As before, 
TT denotes the projection tt : Imm — )■ Bi . Fixing /o and /i , one has 

distf;(7r(/o),7r(/i)) = inf ||/o o - , 



where the infimum is taken over all if e Diff(Af). The Frechet distance is related 
to the Finsler metric 

: TImm(Af,M") ^ M, ^^||^^|Loo- 

Lemma. The pathlength distance induced by the Finsler metric G°° provides an 
upper hound for the Frechet distance: 

dist^{Fo,Fi) < distal. {Fo,Fi) = inf / \\ft\\a^ dt, 

f Jo 

where the infimum is taken over all paths 

/: [0,l]^Imm(Af,E") with 7r(/(0)) = Pq, ^(/(l)) = J^i- 

Proof. Since any path / can be reparametrized such that ft is normal to /, one has 



/Vt^lLoo dt = mi f'wfth^dt, 
J Jo J Jo 



where the infimum is taken over the same class of paths / as described above. 
Therefore 



dist^-(Po, F^) = inf 11/(1) - /(0)||^^ = inf 



/ ftdt <inf f \\fth..dt 

Jo L=° ' JO 

iiif / ||/,^||^^d< = dist^^(Fo,Fi). □ 
/ Jo 
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It is claimed in 15^ theorem 13] that doo — distcoo . However, the proof given 
there only works on the vector space C°°{M, M") and not on Bi{M, M"). The reason 
is that convex combinations of immersions are used in the proof, but that the space 
of immersions is not convex. 

9.4. Sobolev versus Prechet distance. It is a desirable property of any distance 
on shape space to be stronger than the Frechet distance. Otherwise, singular points 
of a shape could move arbitrarily far away without increasing the distance much. 

As the following result shows, Sobolev metrics of low order do not have this 
property. The authors believe that this property is true when the order of the 
metric is high enough, but were not able to prove this. 

Lemma. Let be a metric on _Bj(A/, M") that is weaker than or at least as weak 
as a Sobolev -metric with p < + 1, i.e. 



Then the Frechet distance can not be bounded by the -distance. 

Proof. It is sufficient to prove the claim for P = l + A^*. Let /o be a fixed immersion 
of AI into M", and let /i be a translation of /o by a vector h of length £. It will be 
shown that the iJ^-distance between 7r(/o) and 7r(/i) is bounded by a constant 2L 
that does not depend on £, where tt denotes the projection of Imm onto Bi. Then 
it follows that the i/^'-distance can not be bounded from below by the Frechet 
distance, and this proves the claim. 

For small rg, one calculates the _ff^'-length of the following path of immersions: 
First scale /o by a factor tq, then translate it by h, and then scale it again until 
it has reached /i. The following calculation shows that under the assumption 
p < to/2 + 1 the immersion /o can be scaled down to zero in finite iJ^-pathlength 
L. Let r : [0, 1] [0, 1] be a function of time with r(0) = 1 and r(l) = 0. 



The last integral converges if ^ < which holds by assumption. Scaling 
down to rg > needs even less effort. So one sees that the length of the shrinking 
and growing part of the path is bounded by 2L. 

The length of the translation is simply ^\/''™Vol(/o) = 0(r™/^) since the Lapla- 
cian of the constant vector field vanishes. Therefore 




for all h G Tlmm. 





□ 
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9.5. Concentric spheres. For a Sobolev type metric that is invariant under 
the action of the SO{n) on M", the set of hyper-spheres in M" with common center 
is a totahy geodesic subspace of 5^(5'"^^, M"). The reason is that it is the fixed 
point set of the group SO{n) acting on Bi isometricahy. (One also needs uniqueness 
of solutions to the geodesic equation to prove that the concentric spheres are totally 
geodesic.) This section mainly deals with the case P = 1 + A^. 

First we want to determine under what conditions the set of concentric spheres 
is geodesically complete under the G^-metric. 

Lemma. The space of concentric spheres is complete with respect to the metric 
with P = 1 + AAP iffp>in+ l)/2. 



Proof. The space is complete if and only if it is impossible to scale a sphere down 
to zero or up to infinity in finite G^ path-length. So let / be a path of concentric 
spheres. It is uniquely described by its radius r. Its velocity is ft = rt.v, where v 
designates the unit normal vector field. One has 



g{g-\S,,^)=:L = -^ldTM, Tr(L^) = (-1)'=^^, Yol = r^-'^ 



r(n/2+l)- 

Keep in mind that r and r* are constant functions on the sphere, so that all deriva- 
tives of them vanish. Therefore 



Aiy = V*(V//) = V*{-Tf.L) = Tr^ [V{Tf.L) 
= Tr^ (v(T/).i) + Trf (t/.(VL)) 
= Tt{L^).i. + Tr9 fr/.V(-i Wtm)) - ^.^^ + 



Pft^{l + AAP){rt.,^)^rt.[l + A^ > 



and 



r2p 

From this it is clear that the path / is horizontal. Therefore its length as a path in 
Bi is the same as its length as a path in Imm. One calculates its length as in the 
proof of|9.4| 



(n- l)P 



vo\{g)dt 



{n~l)P\ n.7r"/2 



V2 / (n-l)f 



y r(n/2 + 1) y V r^p 

The integral diverges for ri — > oo since the integrand is greater than r^"^^-'/^. It 
diverges for tq — ;> iff (n— 1 — 2p)/2 < —1, which is equivalent top > (n + l)/2. □ 
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The geodesic equation within the space of concentric spheres reduces to an ODE 
for the radius that can be read off the geodesic equation in section [873) 

2/71-1 p.A.{n-l)P 

~ V 2r " r{r^P + A{n-l)P) 

10. DiFFEOMORPHISM GROUPS 

For M = N the space Emb(M, Af) equals the dijjeomorphism group of M. An 
operator P S r(L(TEmb; TEmb)) that is invariant under reparametrizations in- 
duces a right-invariant Riemannian metric on this space. Thus one gets the ge- 
odesic equation for right-invariant Sobolev metrics on diffeomorphism groups and 
well-posedness of this equation. To the authors knowledge, well-posedness has so 
far only been shown for the special case M = 5^ in [7] and for the special case of 



Sobolev order one metrics in [lOj . Theorem 6.6 establishes this result for arbitrary 



compact M and Sobolev metrics of arbitrary order. 

In the decomposition of a vector h £ T^Emb into its tangential and normal 
components h — Tf.h^ + h-^, the normal part h-^ vanishes. Also S = VT/ vanishes. 
Thus the geodesic equation on DifF(M) in terms of the momentum p is given by 



(see 6.5) 



j p = Pft(^yo\{g), 
\ Vg^p = -Tf.g{PfuV0 ® vol(5). 
Note that this equation is not right-trivialized, in contrast to the equation given in 



[Tl[T8l[T6], for example. The special case of theorem 6.6 now reads as follows: 
Theorem. Let p > 1 and k > '^'™(-'''^) ^ i djid l^i p satisfy assumptions (1-3) of 

Then the initial value problem for the geodesic equation has unique local solutions 
in the Sobolev manifold DiS'^^^^ of H^'^^p -diffeomorphisms. The solutions depend 
smoothly on t and on the initial conditions /(O, . ) and ft(0, . ). The domain of 
existence (in t) is uniform in k and thus this also holds in Diff (Af ) . 

Moreover, in each Sobolev completion Diff'^^'^^, the Riemannian exponential 
mapping exp^ exists and is smooth on a neighborhood of the zero section in the 
tangent bundle, and (7r,exp^) is a diffeomorphism from a (smaller) neigbourhood 
of the zero section to a neighborhood of the diagonal in Diff*^^^^ x Diff'^^^'". All 
these neighborhoods are uniform in k > dim(Af)/2 -I- 1 and can be chosen H^o+'^v. 
open, for fcg > dim(Af)/2 -I- 1. Thus both properties of the exponential mapping 
continue to hold in Diff(Af). 



11. Numerical results 



It is of great interest for shape comparison to solve the boundary value problem 
for geodesies in shape space. When the boundary value problem can be solved, 
then any shape can be encoded as the initial momentum of a geodesic starting at 
a fixed reference shape. Since the initial momenta all lie in the same vector space, 
this also opens the way to statistics on shape space. 
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There are two approaches to solving the boundary value problem. In [2] the 
first approach of minimizing horizontal path energy over the set of curves in Imm 
connecting two fixed boundary shapes has been taken. This has been done for 
several almost local metrics. For these metrics it is straightforward to calculate 
the horizontal energy because the horizontal bundle equals the normal bundle. 
However, in the case of Sobolev type metrics the horizontal energy involves the 
inverse of a differential operator (see section 6.8 1, which makes this approach much 
harder. 





Figure 1. Geodesic where a bump is formed out a flat plane. 
The initial momentum is a = sin(a;) sin(y). Time increases linearly 
from left to right. The final time is t = 5. 



0.8250 h 
0.8245 
0.8240 
0.8235 
0.8230 



0.8225 




along the geo- 
0.822467 



Figure 2. Conservation of the energy G^{ft,ft 
desic in figurejlj The true value of G^{ft, ft) is 4(^i\2A) 
for A ~ 1. The maximum time step used in blue and green is 0.1. 
For purple and cyan it is 0.05. The number of grid points used in 
blue and cyan is 100 times 100. For green and purple it is 200 
times 200. 



The second approach is the method of geodesic shooting. This method is based 
on iteratively solving the initial value problem while suitably adapting the initial 
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-^^^^ -^^^^ ^^^^ 



Figure 3. Letter A forming along a geodesic path. Time increases 
linearly from left to right. The final time is < = 0.8. Top and 
bottom row arc different views of the same geodesic. 




Figure 4. Initial velocity /t(0, •) and momentum a(0, •) of the 
geodesic in figure [3j Both are shown first from above, then from 
the side. 




Figure 5. A self-intersection forming along a geodesic. Time in- 
creases linearly from left to right. 

conditions. The theoretical requirements of existence of solutions to the geodesic 
equation and smooth dependence on initial conditions are met for Sobolev type 
metrics, see section [6. 6| This makes geodesic shooting a promising approach. 

The first step towards this aim is to numerically solve the initial value problem 
for geodesies, at least for the -metric and the case of surfaces in M."^ , and that is 
what will he presented in this work. 

The geodesic equation on shape space is equivalent to the horizontal geodesic 
equation on the space of immersions. For the case of surfaces in E^, it takes the 
form given in section |9.1| This equation can be conveniently set up using the 
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DifFerentialGeometry package incorporated in the computer algebra system Maple 
as demonstrated in figure |6] (The equations that have actually been solved were 
simplified by multiplying intermediate terms with suitable powers of y^Yol{g), but 
for the sake of clearness this has not been included in the Maple code in figure [6]) 



1 with(Dif f erentialGeometry) ; with (Tensor) ;with(Tools) ; 

DGsetupC [u, v] , [x,y ,z] ,E) ; 
3 f := evalDG(f l(t,u,v)*D_x+f2(t,u,v)*D_y+f3(t,u,v)*D_z) ; 

G := evalDG(dx &t dx + dy &t dy + dz &t dz) ; 
5 Gaimna_vrt := femult Connection(dx &t D_x &t du) ; 

Tf := CovariantDerivative(f ,Gainma_vrt) ; 
7 g:=ContractIndices(G &t Tf &t Tf , [ [1 ,3] , [2 ,5] ] ) ; 

g_inv : =InverseMetr ic (g) ; 
9 Gaimna_bas : =Christof f el (g) ; 

det : =Hook ( [D_u , D_u] , g) *Hook ( [D_ v , D_v] , g) -Hook ( [D_u , D_v] , g) " 2 ; 
11 ft := evalDG(diff (f l(t,u,v) ,t)*D_x+diff (f2(t,u,v) ,t)*D_y 

+diff (f3(t,u,v) ,t)*D_z) ; 
13 ft := convert (f t, DGtensor) ; 

S := CovariantDerivative (Tf ,Gainma_vrt ,Gamma_bas) ; 
15 cross :=evalDG((dy few dz) &t D_x + (dz &w dx) &t D_y + (dx &w dy) &t D_z) ; 

N:=Hook( [ContractIndices(Tf &t D_u, [ [2 ,3] ] ) , 
17 ContractIndices(Tf &t D_v, [ [2,3] ])], cross) ; 

nu:=convert(N/sqrt(Hook( [N,N] ,G)) .DGvector) ; 
19 s := ContractIndices(G &t S &t nu, [ [1 ,3] , [2,6] ] ) ; 

L := Contractlndices (g_inv &t s,[[2,3]]); 
21 Gftft := Contractlndices (G &t ft &t f t , [ [1 , 3] , [2 ,4] ] ) ; 

Cft := CovariantDerivative (ft ,Gamina_vrt) ; 
23 CCft := CovariantDerivative (Cf t, Gaiimia_vrt, Gaiimia_bas) ; 

Dft := Contractlndices (-g_inv &t CCf t , [ [1 ,4] , [2, 5] ] ) ; 
25 Pft := evalDG(ft+A*Dft) ; 

GCftCft := Contractlndices (Cft &t Cft &t G, [ [1 ,5] , [3, 6] ] ) ; 
27 gGCftCft := Contractlndices (g_inv &t GCf tCf t , [ [2 ,3] ] ) ; 

TrLgGCftCft := Contractlndices (L &t gGCf tCf t , [ [1 ,4] , [2,3] ] ) ; 
29 TrgGCftCft := Contractlndices (gGCf tCft ,[ [1 , 2] ]) ; 

TrL := Contractlndices (L, [ [1 , 2] ]) ; 
31 # b(t,u,v) := a(t,u,v)*sqrt(det) ; 

eql := Contractlndices (Pft &t dx, [[1,2]] )*sqrt(det) = 
33 b(t,u,v)*ContractIndices(nu &t dx,[[l,2]]); 

eq2 := Contractlndices (Pft &t dy , [ [1 ,2] ] )*sqrt (det) = 
35 b(t,u,v)*ContractIndices(nu &t dy,[[l,2]]); 

eq3 := Contractlndices (Pft &t dz , [ [1 ,2] ] ) *sqrt (det) = 
37 b(t,u,v)*ContractIndices(nu &t dz,[[l,2]]); 

eq4 := diff (b(t ,u, v) ,t) = 
39 (A*TrLgGCftCft - TrL/2* (Gf tf t+A*TrgGCf tCf t) ) *sqrt (det) ; 



Figure 6. Maple source code to set up the geodesic equation. 
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Unfortunately, Maple (as of version 14) is not able to solve PDEs with more than 
one space variable numerically. Thus the equations were translated into Mathemat- 
ica. The PDE was solved using the method of lines. Spatial discretization was done 
using an equidistant grid, and spatial derivatives were replaced by finite differences. 
The time-derivative ft appears implicitly in the equation Pfift) = a-*^, and this 
remains so when the operator Pf is replaced by finite differences. 

The solver that has been used is the Implicit Differential- Algebraic (IDA) solver 
that is part of the SUNDIALS suite and is integrated into Mathematica. IDA 
uses backward differentiation of order 1 to 5 with variable step widths. Order 5 is 
standard and has also been used here. At each time step, the new value of ft is 
computed using some previous values of /, and then the new value of / is calculated 
from the equation Pfift) — a.v. The dependence on / in this equation is of course 
highly nonlinear. A Newton method is used to solve it. This operation is quite 
costly and has to be done at every step, which is a main disadvantage of backward 
differentiation algorithms. Explicit methods are probably much better adapted to 
the problem. The implementation of an explicit solver is ongoing common work of 
the authors with Martins Bruveris and Colin Cotter. 

In the examples that follow, / at time zero is a square [0, tt] x [0, tt] flatly em- 
bedded in M.^ . This is a manifold with boundary, but it can be seen as a part of 
a bigger closed manifold. Zero boundary conditions are used for both / and a. 
It remains to specify an initial condition for a. As a first example, let us assume 
that a at time zero equals sin(a:) sin(?/), where x^y are the Euclidean coordinates 
on the square. The resulting geodesic is depicted in figure [T] In the absence of a 
closed-form solution of the geodesic equation, one way to check if the solution is 
correct is to see if the energy G^{ft,ft) is conserved. Figure [2] shows this for the 
geodesic from figure [l] with various space and time discretizations. 

A more complicated example of a geodesic is shown in figure [3]and[4l There, the 
initial velocity was chosen to be a smoothened version of a black and white picture 
of the letter A. The initial momentum was computed from it using a discrete Fourier 
transform. 

Finally, it is shown in figure [5] that self-intersections of the surface can occur. 
This is not due to a numerical error but part of the theory, and can be an advantage 
or a disadvantage depending on the application. 
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